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Modifications made in the thesis titled 

“BAND-LIMITED WAVELETS AND WAVELET PACKETS” 


by 

Bis^aranjan Behera 

1. Page 6, lines 5-6. The sentence 

“In fact, such wavelets exist in plenty and are easy to construct”, 
should be changed to the following: 

“Such wavelets exist in plenty and are easy to construct. In fact, the 
Shannon wavelet, presented in Example 1.2, gives an affirmative answer 
to the question”. 

2. Page 8. Corollary 1.2 should be rewritten as follows: 

Corollary 1.2. Let K and W be subsets ofM., and W is both transla- 
tion and dilation equivalent to K . Then W is a wavelet set if and only 
if K is so. 

3. Page 10. The last sentence of section 1.3 should be changed to: 

“In Chapter 2, we construct a large class of non-MRA wavelets; the 
two examples w r e just presented belong to this class” . 

4. Page 16. Remark 2.1.1 is incorrect and should be, therefore, deleted. 

o. Page 24. The first sentence should be changed to the following: 

“We now present a proof of the last statement in Theorem 2.1”. 

6. At the end of Chapter 4, the following remark should be added. 

Remark 4.1. G. Majchrowska [Maj] has also constructed some new 
wavelet sets of J3*(R). In particular, she proves the existence of 
an if 2 -wavelet which does not belong to any U J (W), p < 2. 

I am grateful to the thesis examiners for pointing out the reference [Maj] as 
well as for the suggestions for the above modifications in the thesis. 

Additional Reference 

[Maj] Majchrowska, G. (2001). Some new examples of wavelets in the Hardy 
space H 2 (R) . Bulletin of the Polish Academy of Sciences: Mathematics, 
Vol. 49, No. 2. 
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Synopsis 

In this thesis we construct large classes of band-limited wavelets of L 2 (R) and the Hardy 
space H 2 (R). We have also considered the construction of multiwavelet packets and 
wavelet frame packets of L 2 (R d ) with arbitrary dilation matrices. 

The Fourier transform of a function / € X 1 (R <i ), d > 1 is defined by 

mo = /({) = [ f 6 I s . 

Ju d 

To define the Fourier transform for functions of L 2 ( R d ), the operator T is extended from 
L 1 n L 2 (B. d ) to the whole of L 2 (R d ). A function / is said to be band-limited if / is 
compactly supported. 

Definition 0.1. A function ip E L 2 { R) is said to be a wavelet of L 2 (R) if the system of 
functions { ip jt k : j,k €%} forms an orthonormal basis for L 2 (R), where 

iPj, k (x) = 2^(yx-k), j, k£l. 

G. Gripenberg and X. Wang independently proved the following theorem which char- 
acterizes all wavelets of L 2 (R) . 

Theorem 0.1. A function ip, with ||^>||2 = 1 , is a wavelet for L 2 (R) if and only if 

53 IV>( 2i £)l 2 = 1 for a - e - £ € R 

j€Z 

53 $$£)'$ (2 J (£ + 2 ? 7 r)) = 0 for a.e. £ € R and for all q € 2Z + 1. 

j>0 

Y. Meyer and S. Mallat defined the notion of a multiresolution analysis, which provides 
a framework to construct wavelets of L 2 (R). It is defined as follows. 
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Definition 0.2. A multiresolution analysis (MRA) is a sequence of closed subspaces 
[V 0 : j G Z} of I? (R) satisfying the following properties: 

(i) V 0 C Vj+i for all j G Z 

(ii) UjzzVj is dense in L 2 ( R) and fljgzV^ = {0} 
fmj / eVj if and only if /( 2-) G V^+i /or all j G Z 

(iv) there exists a function G L 2 (R) suc/i i/iai {(/?(• — k) : A: G Z forms an orthonormal 
basis for Vq. 

Given a wavelet if, define the subspaces V 3 = sp{ipi,k '• l < j, k G Z}. If {1^ : j € Z} 
is an MRA, then we say that ip is an MRA-wavelet. Let 

]>i kez 

The function D ^ is called the dimension function of the wavelet ip. The following char- 
acterization of MRA-wavelets was proved independently by Gripenberg and Wang. 


Theorem 0.2. A wavelet ip is an MRA wavelet if and only if D^, = 1 a.e. 


Meyer constructed a wavelet such that both the wavelet and its Fourier transform 
are infinitely differentiable. The support of the Fourier transform of Meyer s wavelet is 
[-87,. _?j u [lyr I;,-], e. Hernandez, X. Wang, and G. Weiss characterized all wavelets 
with Fourier transform supported in this set (see Theorem 4.1, Chapter 3 of [HW]). Our 
first result is the following. In Chapter 2, for n > 1, we construct a set S n ? which is 
symmetric with respect to the origin. The set Si is the Shannon set [—2 7r, — 7 t] U [tt, 27t] 
whereas S 2 is the set [— fvr, — |tt] U [|tt, §tt] associated with the Meyer wavelet. Then we 
characterize all wavelets with Fourier transform supported in S n . Let n > 1. Put 


2^1— i 2 n 

° n = bn = 2^1 


2 n ~ 1 ( 2 n — 2 ) 2 2n ~ 1 

tt, Cn = — — — ; 7r, d n = — -ft, and e n 


2 n — 1 


2 n — 1 


2 n — 2 
2 n - 1 


7 r. 


(2) 


Define 5„ = 5+ U where S+ = [a n , b n ] U [c*, d n ] and S n = -(S+). 
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We prove the following theorem which characterizes all wavelets xp such that 
supp ip C S n . 

Theorem 0.3. Let n > 2, ip e L 2 (R) , supp xp C S n and b(£) = |$(£)|. Then ip is a 
wavelet for L 2 (R) if and only if 

(i) &(£) 1 for a.e. £ £ [<2n> e n ] U [ e n , a n ] 

(ii) f> 2 (£) + b 2 ( 2 n ~'Z) = 1 for a.e. [e n , b n ] 

(in) i> 2 (£) 4- b 2 (£ - 2?r) = 1 for a.e. £ € [e n , b n ] 

(iv) 6(0 = b (2 n-1 (f - 2 t r)) for a.e. £ € [e n , b n \ 

(v) xp(£) = e t6 ^b( 0, where 9 satisfies 

0(0 + 9 (2 n - 1 (e - 20) - e{t - 2tt) - 0(2 n-1 O = (2m(0 + 1) ir, 

for some m { 0 £ Z, for a.e. £ £ [e„, 6 n ] n (supp 6) n (j^rsupp &)• 

Moreover, for n > 3, none o/ iftese wavelets arises from an MRA. 

For n > 3, we have also computed the dimension functions D n , defined in (1), of 
wavelets associated with the set S n . These dimension functions are even and attain all 
integral values from 0 to the maximum value attained by D n . 

Chapter 3 deals with the constructions of wavelets with interesting properties. These 
wavelets are constructed with the help of the wavelets associated with the set S n of 
Chapter 2. Our first result in this chapter is the construction of a family of band-limited 
wavelets whose Fourier transforms do not vanish in any neighbourhood of 0. 

For n > 2, let a n , b n , Cn, d n , and e n be as in (2) and let 0 < e < g f. Define 

p i=[j + ^j + '}’ a = [« n + 2e,t„], and ft = [4,4 + 2a]. 

Let 

X 0 = Pi— 2-2"- 2 tt, Yo = ±X 0 , 

X l = Y l - l -2-2 n ~ 2 7c, l > 1. 
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Now define the set 

W n = \ jj Xi j U (0.) UMU (P x u P 2 U P^. (3) 

Theorem 0.4. For n > 2, the function ip n , defined by ip n = Xw n > IS a band-limited 
wavelet such that ip does not vanish in any neighbourhood of 0. 


In [Gar] G. Garrigos asked whether there exists a wavelet if) such that (i) if) is band- 
limited, i.e., ip has compact support, (ii) ip is even, and (iii) ip does not vanish in any 
neighbourhood of the origin. Examples of wavelets satisfying any two of the above three 
properties exist in the literature. The next result is the construction of a family of wavelets 
satisfying all the three properties (i)-(iii). For e > 0 such that 0 < e < construct the 
following sets: 



. ^ . 
+ F’ Y + 


S 2 — a n 2e, 7r , — [d n > d n -I- 2e] . 


Let T t = —S t for i = 1, 2, 3. Now, let 


E 0 = S 1 + 2-2 n 2 7t, Fo = 

Ei = F/_i + 2 • 2 n—2 7r, F, = 2 ^+iE h l > 1, 


G[ = -Ei and Hi = —Fi, for l > 0. 


Define 


00 


OO 


K n = (P 2 - |J Eij U (JJ Fij U (Sx US 2 U Si) 

u (l 2 - u (u«) U (Fi UT 2 U Ts). 


0 


U=0 


(4) 


Theorem 0.5. For each n> 2, the function w n , defined by w n = xk u > a wavelet which 
satisfies the above three properties. 


In [Web] E. Weber defined an equivalence relation on the set of all wavelets of L 2 (R) in 
the following manner. For a wavelet ip, define the subspaces V 3 = sp{ip 3 ,k • l < j, k G Z}. 
It is easy to verify that Vq is invariant under translations by integers. Are there other 
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groups of translations under which V 0 is invariant? For a € R, let T a denote the unitary 
operator T a f = /(• — o>). Consider the group of translations Q n = {T™ : m € Z} and 
the group Goo = {T a : a € R}. Let C n be the collection of all wavelets such that the 
corresponding space Vo is invariant under the group Q n . Then we have 

£ 0 => £ x D £ 2 D • • -C n D C n+ i D • • o £ oo. 


An equivalence relation can now be defined on the collection of all wavelets of L 2 (R). 
The equivalence classes are given by M n = £„ \ £ n+1 , with = £00 • Therefore, M n 
is the class of wavelets such that Vo is invariant under the group Q n but not under Gn+i- 
Weber characterized the classes M. n and constructed examples of wavelets in the classes 
Mt, 2 = 0 , 1,2,3. We have constructed examples of wavelets in each of the equivalence 
classes M n , n > 0 . 

Let n > 3 and a n , b n , c n , d n , and e n be as in ( 2 ). Define ip n by 


1 


&(0 = 


75 

1 

~V2 


0 


if ^ G [ — 671 , —On] u [Cn, ^ + 2 n_1 7r] U [^ + 2 n_1 7T, d n ] 
if f e l 9 ?, Tf-] U [o n , e n ] u [^ + 2"~ 1 7t, ^ + 2 n_1 7r] 
if £ € [a„ + 2 n 7T , e n + 2 n 7r] 
otherwise. 


We prove the following theorem. 


Theorem 0 . 6 . The function ip n , defined above, is a wavelet and belongs to the equivalence 
Class Mn-2- 


This theorem proves that M n , n > 1 are non-empty. We have also constructed a 
family of wavelets in the class Mo- 

Chapter 4 is on the construction of wavelets for the Hardy space H 2 ( R) defined by 

H 2 {R) = {fe L 2 (R) : fa) = 0 for a.e. t < 0}. 

A function if € H 2 (R) is said to be a wavelet of H 2 { R), or an if 2 - wavelet, if the system 
of functions {tpj^ := 2 J / 2 ( 2 J • — k ) : j,k € Z} forms an orthonormal basis for H 2 (E) . The 
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function ip such that 4> = X[2*m ] one such example. In fact, this was the only known 
# 2 -wavelet for a long time. A subset K of R+ is said to be an H 2 -wavelet set if Xk 
is the Fourier transform of an FT 2 -wavelet. An H 2 - wavelet set will be called an interval 
H 2 -wavelet set if it is the union of a finite number of intervals. In Chapter 4, we prove 
a result on the structure of interval Ff 2 -wavelet sets and characterize all if 2 -wavelet sets 
consisting of three disjoint intervals. 

In Chapter 5, we extend the notion of wavelet packets, introduced by Coifman, Meyer 
and Wickerhauser, to L 2 (R d ) with arbitrary dilation matrices and we have considered the 
case of the MRAs with multiplicity, i.e., for which the resolution spaces are spanned by 
more than one scaling function. The wavelet frame packets of [Che] is also generalized to 
L 2 (R d ) with arbitrary dilations. 

Let A be a dilation matrix , i.e., A is a d x d matrix such that A(Z d ) C Z d , and all 
eigenvalues of A have absolute value greater than 1. For such a matrix A, | det A\ is an 
integer greater than 1. Let B = A 1 , the transpose of A and a = | det A\ = | det£|. It is 
well known that the number of distinct cosets of AZ d in Z d is a. A subset K of Z d which 
consists of exactly one element from each of these cosets will be called a set of digits for 
A. 

Using the properties of MRA, we can find matrices if r (£), 0 < r < a — 1 of order L, 
having 27rZ d -periodic functions as entries, such that the following relation holds: 

J2 #r(£ + 2£~ V)# s *(f + 2 £" V) = S rs I L for a.e. f, 0 < r, s < a - 1, 

fjtGK 

e I L is the identity matrix of order L, M*(Q is the conjugate transpose of the matrix 
), and K is a set of digits for B. Let 

= 0 < r < a — 1, 

hU() = £ where h T ljt € €. For 0 < r < a - 1, 1 < 1 < L, define 

k€Z d 

f[(x) = E E tfjk al/2( Pj(Ax - k )• 

j = 1 fc€ Z d 
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Now, for any n € No = N U {0}, we define /", 1 < l < L recursively as follows. Suppose 
that f[, r 6 No, 1 <1 < L are defined already. Then define 

ir r (x) = EE h s hk a}l 2 f](Ax - k)\ 0 < s < a - 1, 1 < l < L. (5) 

j=i kez d 

Note that (5) defines /" for every non-negative integer n and every l such that 1 < l < L. 
The functions {/" : n > 0, 1 < l < L} will be called the basic multiwavelet packets 
corresponding to the MRA {Vj}. The first result regarding the wavelet packets is the 
following. 

Theorem 0.7. Let {/" : n > 0, 1 < l < L} be the basic multiwavelet packets. Then 

(i) {/"(• — k) : a 3 < n < a J+1 — 1, 1 < l < L, k € Z d } is an orthonormal basis of 
W v j> 0. 

(ii) {/”(• — k) : 0 < n < a 3 — 1, 1 < l < L, k € Z d } is an orthonormal basis of 

V„ j > 0. 

(iii) {/”(• — k) : n > 0, 1 < l < L, k € Z d } is an orthonormal basis of L 2 (R d ). 

Next, we have considered all dilations by the matrix A and all Z d - translations of the 
basic multiwavelet packet functions. 

Definition 0.3. Let {/” : n > 0, 1 < l < L} be the basic multiwavelet packets. The 
collection of functions V — {a J / 2 /"(A J • —k) : n > 0, 1 < l < L, j G Z, k € Z d } will be 
called the general multiwavelet packets associated with the MRA {Vj} of multiplicity L. 

In the following theorem we characterize all subcollections of V which give rise to 
orthonormal bases of L 2 (R d ). Let No = N U {0} and I no — {r € No : a 3 n < r < 
a 3 (n + 1) — 1}. 

Theorem 0.8. Let {ff : n > 0, 1 < l < L} be the basic multiwavelet packets and 
S C No x Z. Then the collection of functions {atf^A 3 — k) : 1 < l < L, k € Z d , ( n,j ) G 5} 
is an orthonormal basis of L 2 (R d ) if and only if {I n ,j '• i n ij) € 5} is a partition of Nq . 
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A sequence {xk : k € Z} of a separable Hilbert space H is said to be a frame for H if 
there exist constants C\ and C2, 0 < C\ < C 2 < 00 such that for all x e % 

CiW 2 <El^.^)l 2 <c 2 M 2 . 

ke z 

Similar to the orthogonal case, wavelet frame packets can be constructed. This was done 
for L 2 ( E) in [Che] and for L 2 (R d ) in [LC] for the dilation 2. We have extended these 
results to L 2 (R d ) with arbitrary dilations and proved results similar to Theorem 0.7 and 
Theorem 0.8 in this setting. 
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Chapter 1 


Introduction 

1.1 A brief history 

A function is said to be a wavelet if the collection of functions obtained by dyadic dilation 
and integral translations forms an orthonormal basis for L 2 (R). The first example of such 
a function was introduced by Haar in 1909. The disadvantage with this wavelet is that 
it is discontinuous. In 1927 Ph. Franklin constructed a system of continuous orthogonal 
functions which formed a basis of L 2 [ 0, 1]. J.O. Stromberg [Str] introduced the modified 
Franklin system in 1980. The Stromberg wavelet has the following properties: (i) it is 
continuous on the whole real line R; (ii) it is linear on the intervals [0, |], [|, 1] and on 
each interval of the form [k, k + l],k > 1 and of the form [— — |], l > 0; and (iii) 
ip decreases rapidly at infinity. Morlet and Grossman used wavelets in the beginning of 
1980’s for their geological research. Y. Meyer [Meyl] constructed an example of a wavelet 
such that both the wavelet and its Fourier transform are infinitely differentiable. This 
was generalized to the higher dimensional spaces L 2 (R d ) by Lemarie and Meyer in [LM1] 
In 1985 Meyer [Mey2], in collaboration with S. Mallat [Mai], introduced the concept of 
a multiresolution analysis which provided a framework to construct wavelets. Using this 
framework I. Daubechies in 1988 constructed, for each positive integer r, a wavelet which 
is compactly supported and r times continuously differentiable. 


1 
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Let T denote the unit circle in the complex plane which can be identified with the 
interval [— 7T, 7 t). Let L 2 ( T) be the space of all measurable, square-integrable functions 
defined on T. It is a Hilbert space with the inner product 

1 /' 2ir 

(fi9) ~ 2 tt J Q f( x )9(x)dx. 

The set of functions {e m ‘ : n E Z} forms an orthonormal basis for L 2 (T) . By denoting 
g(x) = e ix and g n (x) = g(nx), we see that L 2 (T) is generated by the integral dilations 
of a single function g. That is, any function / in L 2 (T) can be represented as an infinite 
linear combination of the functions g n : 

/(*) ='%2 c n(f)9n(x), 

where the series on the right hand side converges in the norm of L 2 (T), which means that 

n 2 21 r 2 * | 2 

/- X 2 := 2 ^/ I f( x )- X dx — >-0 as N U N 2 — * oo. 

n=—Ni n=—Ni 

The constants Cn(/) are the Fourier coefficients of / defined by Cn(f) = ( f,g n )■ 

If we want to analyse an arbitrary function / of L 2 ( R), the space of square-integrable 
functions defined on the real line R, with the help of functions generated by a single 
function, then such a function has to decay to 0 at ±oo (obviously the functions g n ’s cannot 
be used, as these functions themselves do not belong to L 2 (R)). Therefore, we are looking 
for a small wave, or a wavelet. Since this wavelet will essentially be concentrated on a 
finite interval, to analyse an arbitrary function of L 2 (R) we have to shift the wavelet along 
the real line. That is, the translations of the wavelet have to be considered. Moreover, 
if a function oscillates too often (or, is spread out) in an interval which is very small 
(respectively, large) in comparison to the interval in which the wavelet is concentrated, 
then in order to analyse this function suitably, one has to contract (respectively, stretch) 
the wavelet upto a suitable order. That is, the wavelet has to be dilated. In other words, 
we are looking for a function ip, which is smooth, localized and is such that the set of the 
following dilations and translations 


{ip(2 3 • - k ) : j, k € Z} 
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suffices to analyse all functions of L 2 ( R), as in the case of I» 2 (T). 

1.2 Basic concepts 

In this section we list some of the definitions, examples and results which are well known 
in the theory of wavelets. We start with the precise definition of a wavelet. 

Definition 1.1. A function if G L 2 (R) is said to be an orthonormal wavelet if the system 
of functions { ip h k : j,k G Z } forms an orthonormal basis for L 2 (E), where 

ifj,k(x) = 2 3/2 ip{2 :i x - k), j,kG Z. (1.1) 

In other words, the functions {ip^k} are orthonormal: 

(V’j.fc) = dj,j'bk,k' ) 

and any function / G L 2 (R) can be written as 

jez ke z 

where the convergence of the series is in the L 2 (R)-norm. The inner product of two 
functions g and h of L 2 (M) is given by 

(9,h) = / g(x)h(x)dx. 

J R 

Since we will be concerned only with orthonormal wavelets, we will drop the word or- 
thonormal in order to avoid repetitions. So, whenever we say that ip is a wavelet, we 
mean that it is an orthonormal wavelet in the sense of Definition 1.1. The Fourier trans- 
form of a function plays a crucial role in the wavelet theory, which is defined for a function 
/ G L 1 (R d ), d> 1 as 

Ff(() = /K) = [ ( 6 R*. 

JR d 
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To define the Fourier transform for functions of L 2 (R d ), the operator J- is extended from 
L l fl L 2 (R d ), which is dense in L 2 (R d ) in the X 2 -norm, to the whole of L 2 (R d ). For this 
definition of the Fourier transform, Plancherel theorem takes the form 

(/.»>= /,S€L 2 ( b 4 ). 

Note that 

(^) A (0 = 2-i e - l2 ' Jfc ^(2-^), e € R, j, k 6 Z. (1.2) 

Two simple examples of wavelets are the Haar wavelet and the Shannon wavelet. The 
Haar wavelet is compactly supported in the time domain, i.e., in the x-variable and the 
Shannon wavelet is compactly supported in the frequency domain, i.e., in the ^-variable. 

Example 1.1. The Haar wavelet is defined by ip H {x) = X[o,i/ 2 ](^) - X[i/2,i] ( x )- 

Example 1.2. The Shannon wavelet tp s . This wavelet is defined in terms of the Fourier 
transform: 

i> s = Xs, s= [-27T, -7 r] U [tt, 2tt]. 

Observe that for each integer j, 

ke z. 

Let the inner product in L 2 (S) be given by 

(/, 9) = f f(x)g(x)dx. 

Js 

Since : k € Z j is an orthonormal basis of L 2 (S), the collection of functions 

s : fcezjisan orthonormal basis of X 2 (2 J 5). Observe that {2 j S : j € Z} 
is a pairwise disjoint collection of subsets of R whose union is R. Therefore, by the above 
equality, {-T=(^) A \j,k€ Z} is an orthonormal basis of L 2 (R). By Plancherel theorem, 
: j,k e Z} is an orthonormal basis of L 2 (R). In other words, ip s is a wavelet. □ 

A wavelet ip, in particular, has the property that its integral translates form an or- 
thonormal system. For such functions the following result holds. 
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Proposition 1.1. Let f € L 2 (R). Then {/(• - k): k € Z} is an orthonormal system if 

and only if J2 I /(£ + 2k-ir)\ 2 = 1 for a.e £ 6 R. 

kez 

Proof: Suppose that {/(• — fc) : A; € Z} is an orthonormal system. Then we have 

4,o = </,/(•“*)> 

- 

= 1/(01 

- sjf + 

So the /cth Fourier coefficient of the 27r-periodic function ^ |/(- + 2/7r)| 2 is d> fc;0 . Since the 

J€Z 

function which is 1 a.e. also has 5k,o as its fc-th Fourier coefficient, by uniqueness of the 

Fourier series it follows that £ | /(£ + 2lir)\ 2 = 1 for a.e. £ € R. By reversing the above 

le z 

steps we can prove the converse. □ 

The support of a function g is defined as 

supp g = {x € R : g(x) ^ 0}. 

The Lebesgue measure of a set E of R will be denoted by \E\. The following result gives 
information about the measure of the support of / (Corollary 2.4, Chapter 2 of [HW]). 

Corollary 1.1. Let f € L 2 (R) and {/(• — k): k E Z} is an orthonormal system. Then 
| supp f | > 27T. Equality holds if and only if \f\ = xk for some measurable subset K of R. 

Proof: From Proposition 1.1 it follows that |/(£)| < 1 for a.e. (el Therefore, 

2ir = ll/lll = [ |/(£)| 2 # < / . 1 • <*£ = I supp /|, 

JR J(supp f) 

which proves the first part. Now, suppose that | supp /| = 2i r. If |/(£)| < 1 on a set E of 
positive measure, then 

2tt = f _ |/(£) | 2 d£ < | supp f\E\ + \E\ = \ supp f\ = 2? r. 

y (supp /) 
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Hence, we get a contradiction which proves that |/(£)| = 1 on supp f. Conversely, suppose 
that {/(• - k) : k € Z} is an orthonormal system such that |/| = xk- Then 

\K\ = | supp f\ = \\f\\l = 2t r. 

□ 

A natural question to ask is the following: Does there exists a wavelet such that the 
support of its Fourier transform has the minimal measure 27r? In fact, such wavelets exist 
in plenty and are easy to construct. Such wavelets are called minimally supported 
frequency (MSF) wavelets. In [HKLS] these were called unimodular wavelets. 

Definition 1.2. A wavelet ip is said to be an MSF wavelet if l^l = xk for some measurable 
subset K of R. The associated set K is said to be a wavelet set. 

A function ip € L 2 { R) is a wavelet of L 2 (R) if and only if ip satisfies the following four 


conditions: 

y |^(2 J f)| 2 = 1 for a.e. £ € R. (1.3) 

jez 

y V>(2 J £)^ (2-?’(£ + 2qn)) = 0 for a.e. £ G R and for all q G 2Z + 1. (1.4) 

j>0 

y |^(£ + 2&7r)| 2 = 1 for a.e. (£R. (1.5) 

kez 

y ip(£ + 2 kTr)ip (2J(£ + 2 & 7 t)) = 0 for a.e. £ € R and for all j > 1. (1.6) 

fee z 


The fact that equations (1.3)-(1.6) characterize all wavelets of L 2 ( R) was observed by 
many authors. Lemarie and Meyer [LM2] obtained these equations for compactly sup- 
ported wavelets. Bonami, Soria and Weiss [BSW] proved this fact under the assumption 
that ip is band-limited, i.e., ip is compactly supported. For a function ip to be a wavelet, 
it is necessary that ||^|| 2 = 1. With this assumption on ip, the following theorem was 
proved independendly by G. Gripenberg and X. Wang (see also [HKLS]). 

Theorem 1.1. Let ip G L 2 (R) with \\ 1 pW 2 = 1- Then ip is a wavelet of L 2 ( R) if and only 
if ip satisfies (1.3) and (1.4). 
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For a proof of this theorem we refer to Chapter 7 of [HW]. It was observed in [BSW] 
that the orthonormality of the system {xp Jtk :j,ke 1} is characterized by (1.5) and (1.6), 
and completeness by (1.3) and (1.4) (see also [HW]). A consequence of Theorem 1.1 is 
that equations (1.3) and (1.4), along with the assumption = 1, imply the equations 
(1.5) and (1.6). 

For MSF wavelets the characterization is simplified further, in the sense that only 

(1.3) and (1.5) are needed as observed by several authors (see [HKLS], [HW]). 

Theorem 1.2. Let xp G L 2 (R) be such that \xp\ = xk for some measurable subset K of 
R. Then xp is a wavelet if and only if (1.3) and (1.5) hold. 

Proof: Since every wavelet must satisfy (1.3)-(1.6), we need only to show that if xp is 

as in the hypothesis and satisfies (1.3) and (1.5), then it is a wavelet. We will show 

that (1.3) =>(1-6) and (1.5) =^(1.4). If 2 G supp xp = K, then \xp{2 3 ^)\ = 1. As 

Y |^(2 J £ + 2&7r)| 2 = 1 by (1.5), we obtain that 0(2 :? £ + 2A;7r) = 0 for k ^ 0. In particular, 
k€ z 

xp( 2 J £ + 2 • 2 3 qir) = 0 for every odd integer q which proves that each term of the series in 

(1.4) is zero. Now if £ + 2kit € K, then | xp(£ + 2fc7r)| = 1. Writing £ + 2kn instead of £ in 

(1.3),; we have Y lV’(2 J (£ + 2/i;7r))| 2 = 1. So xp(2 3 (£-+2k / n)) = 0 for all j ^ 0. Hence, each 
jez 

term of the series in (1.6) is zero. □ 

The above theorem for MSF wavelets can be restated in terms of the associated wavelet 
sets. 

Theorem 1.3. Let K C R. Then K is a xuavelet set if and only if 

1. {K + 2kn : k € Z} is an a.e. partition of R. 

2. {2°K : j e Zj is an a.e. partition of R. 

We say that {Aj} is an a.e. partition of X if Ai C X for all l, \ Ai D Ai>\ = 0 for l ^ l', 
and \X \ UjA/| = 0. To construct a wavelet set from another wavelet set, the concepts of 
dilation equivalence and translation equivalence of sets will be very helpful. 
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Definition 1.3. A measurable set A is translation equivalent to a measurable set B if 
there exists a measurable partition {A n } of A and k n G Z such that {£„} = {A n + 2k n ir} 
is a partition of B. Similarly, a measurable set A is dilation equivalent to a measurable set 
B if there exists a measurable partition {A n } of A and j n G Z such that {B n } = (2 3n A n } 
is a partition of B. 

As a consequence of Theorem 1.3 we have 

Corollary 1.2. Suppose K is a wavelet set and W is both translation and dilation equiv- 
alent to K. Then W is also a wavelet set. 


1.3 Multiresolution analysis 

The concept of multiresolution analysis (MRA), formulated by Y. Meyer and S. Mallat, 
is crucial to the theory of wavelets. Given an MRA, one can always construct a wavelet 
of L 2 (R). An MRA is defined as follows. 

Definition 1.4. A multiresolution analysis (MRA) is a sequence of closed subspaces 
{V) : j € Z} of L 2 ( R) satisfying the following properties: 

(i) Vj C Vj+i for all j G Z 

(ii) UjgzV; Is dense in L 2 ( R) and = {0} 

(Hi) f €Vj if and only if /( 2-) G V .,+1 for all j € Z 

(iv) there exists a function <p G L 2 (R) such that {< p(- — k ) : k G Z} forms an orthonormal 
basis of Vo. 

A function ip that satisfies property (iv) is called a scaling function for the MRA. If 
Pi is another scaling function for the same MRA, then <pi = T&, where l is a 27r-periodic 
function in L 2 ( T) such that \l(£)\ = 1 a.e. (see Lemma 2.6, Chapter 2 of [HW]). Note 
that {22ip(2 j • -k) : ft G Z} is an orthonormal basis of V 0 because of (in) and (iv). 
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We will now briefly mention how to get a wavelet from an MRA. Since p G Vo C Vi 
and {2^(2 • — k) : A; G Z} is an orthonormal basis of Vy, there exists {c*, : k € Z} € l 2 ( Z) 
such that 

= Y^ ck2 ^v( 2x ~ k )- 

ke 2 

Taking Fourier transforms of both sides, we get 

0(0 = ^ c k^ e ~ lk h (!) = rn 0 (!) <p (!) , (1.7) 


where 

m 0 (O = 2~^J2 Cke ~ tK - 

k€L 

The function m 0 , which is 27r-periodic and is in L 2 (T), is called the low pass filter asso- 
ciated with the scaling function < p . By Proposition 1.1 we have I0(£ + 2 /c 7 t)| 2 = 1. 

k€Z 

Splitting the sum over Z into sums over even and odd integers and making use of (1.7), 
we obtain 

K(0l 2 + K (£ + *■) I 2 = 1 for a.e £ G T. (1.8) 

Let Wj be the orthogonal complement of V 3 in V) +1 : 

w^v^eVj, je z. 

By properties (i) and (ii) of Definition 1.4 we can decompose L 2 (R) as an orthogonal 
direct sum of the Wfi s: 

I 2 (E)=®^. (1.9) 

jez 

It follows from (Hi) that if {ip(- — k) : k G Z} is an orthonormal basis of Wo, then 
{2^(2* • —k) : k e Z} is an orthonormal basis of W y Then by the decomposition (1.9) it 
follows that ip is a, wavelet. Therefore, all that is needed is to find a function ip € Wo such 
that { ip( • - k) : k G Z} is an orthonormal basis of W 0 . Functions ip with this property 
are characterized by 


$(0 = e l 'MO™o (! + tt)0 (!) , 


for a.e (GR, 


( 1 . 10 ) 
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where v is a 27r-periodic measurable function such that |z/(£)| = 1 for a.e. f e T, ip is 
the scaling function of the MRA and mo is the low pass filter associated with <p (see 
Proposition 2.13, Chapter 2 of [HW] for details). 

Given a wavelet ip, define the spaces V 3 = sp{ipi !k :l < j,k € Z}. If {V 3 : j e Z} is an 
MRA, then we say that ip is associated with an MRA (or, ip is an MRA- wavelet) . Not all 
wavelets are associated with an MRA. The first example of a non-MRA wavelet, which is 
an MSF wavelet, was given by J.L. Journe. The associated wavelet set is 


[32 1 


4 ' 


'4 " 


■ 32 

— zr7T, —47 r 
7 

u 

-7T,--7r 

u 

-7T,7r 

u 

4tt, — 7T 

i 


1 


1 


1 


Another interesting example is the Lemarie wavelet which is also an MSF wavelet with 
the following wavelet set: 


'8 4 ' 


‘4 6 ' 


'24 32 ' 


U 

— 7T, — 7T 

U 

— — 7T, —7 r 

7 ’ 7 J 


[7 ’7 J 


7 7 


In Chapter 2 we will show that these two wavelets cannot be constructed from any MRA. 


1.4 The dimension function 


There is a single equation which characterizes all wavelets associated with MRAs. Before 
stating this result let us see how the Fourier transform <p of the scaling function of an 
MRA is related with the Fourier transform ip of the corresponding wavelet. From (1.7), 
(1.8) and (1.10) we get 

|£(2f)| 2 + |^(2£)| 2 = (K(£)| 2 + |m 0 (f + 7r)| 2 ) • |£(£)| 2 = |£(£)| 2 - 


Iterating this equation, we obtain 

lv(?)| 2 = Isp( 2"{)| 2 + E l^ 2 '«)l 2 for aU N 2 1- 


N 


i= i 


Since \<p{£)\ < 1, the sequence j J] \ip(2 j €)\ 2 : N > lj, which is an increasing sequence 


N 


of real numbers, is bounded above by 1 so that lim \ip{2 3 ^)\ 2 exists. Therefore 


N-¥ OO j=1 
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lim |<£(2^f)| 2 also 

N-+00 


exists. Now by Fatou’s lemma, we have 

< Km 


lim -T7 

N-t oo 2 ^ 


27T 


= 0 . 


Hence, 

10(0 1 2 = X) I0( 2J O| 2 for a.e. (1.11) 

j>i 

Using Proposition 1.1, we obtain 

D A0 = £ E 10 ( 2 '(£ + 2 M) I 2 = 1 for a.e. £ € R (1.12) 

j>i fcez 

Therefore, for a wavelet ip to be associated with an MRA, it is necessary that the func- 
tion D^p is equal to 1 a.e. The fact that it is also sufficient is proved independently by 
Gripenberg and Wang. We refer to Chapter 7 of [HW] for a proof of the following theorem. 

Theorem 1.4. A wavelet ip is an MRA wavelet if and only if = 1 for a.e. £ € R. 

The function D ^ is called the dimension function of ip, since for a.e. £, A/,(£) happens 
to be the dimension of certain subspaces of l 2 ( Z) (see [Aus]). This fact, in particular, shows 
that (£) is an integer for a.e. £. The significance of the dimension function is that, 
given a dimension function one can construct wavelets associated with it. In [BRS], a 
characterization of dimension functions is given. An algorithm to construct wavelets from 
the dimension function is also provided therein. 


1.5 Wavelets for H 2 (R) 

The space H 2 { R) is the collection of all functions in L 2 (R) whose Fourier transform vanish 
on the negative real axis: 


H 2 { R) = {/ € L 2 (R) : /(£) = 0 for a.e. £ < 0}. 
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It is clear that H 2 (R) is a closed subspace of L 2 (R). As in the case of L 2 (M), we can 
define wavelets for H 2 ( R). A function ip € H 2 (R) is said to be a wavelet for H 2 (R) (or, 
an H 2 -wavelet) if the system of functions {ip 3 ,k : j, k € Z} forms an orthonormal basis for 
H 2 (R). 

Example 1.3. The function ip given by ip = X[2*M is an # 2 -wavelet. The proof of 
this fact is easy and is, in fact, similar to the case of Shannon wavelet of L 2 (R) (see 
Example 1.2). 

Analogous to the L 2 (R) case, we call a set if C (0, oo) an if 2 -wavelet set if xk is 
the Fourier transform of an H 2 - wavelet. If if is a union of a finite number of intervals 
then K will be called an interval wavelet set. In [HKLS] the authors characterized all 
H 2 - wavelets ip such that \ip\ is the characteristic function of either an interval or is the 
union of two disjoint intervals. In Chapter 4 we will prove a result on the structure of 
the interval wavelet sets of H 2 { R) and characterize all such sets which are unions of three 
disjoint intervals. 

1.6 Multiwavelet packets of L 2 (R d ) 

Let us recall that a function ip is a wavelet of L 2 (R) if {2 3 ^ 2 ip(2 3 ■ — k ) : j,k € Z} forms an 
orthonormal basis for L 2 (R). In order to extend this concept to higher dimensional spaces 
L 2 (R d ), d > 1, the notions of dilation and translation have to be generalized. The natural 
generalization is dilation by powers of 2 in every coordinate, i.e., x — (xi,x 2 , . . . ,xf) t-» 
2 3 x = (2 J xi, 2 3 x 2 , . . . , 2 3 xf), and translation by elements of Z d . Then one can define MRA 
for L 2 (R d ). But we can do more: we can replace the dyadic dilation by a general dilation 
matrix (i.e., a d x d matrix which leaves Z d invariant, and such that all its eigenvalues 
have absolute value greater than 1). Moreover, in the definition of MRA, we can allow 
the central space Vo to be generated by more than one scaling function. Such an MRA 
is called an MRA with multiplicity and the wavelets associated with the MRA are called 
multiwavelets. 
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Let {V, : j € Z} be an MRA of L 2 (R) and W 3 = V 0+ i © V, be the corresponding 
wavelet subspaces so that the orthogonal decomposition L 2 (R) = ® 3& W 3 holds. If xp is a 
wavelet associated with this MRA, then {xpj t k = 2^ 2 xp(2 :> • — k ) : k € Z} is an orthonormal 
basis of Wj and consequently, {xp Jt k • j,k € Z} is an orthonormal basis of L 2 (M) . If %p 
is band-limited, i.e., xp has compact support, then \supp xp 3jk \ = 2 J \supp xp | which shows 
that the wavelets have poor frequency localizations (the supports of xp 3t k gets larger for 
large j). In order to overcome this disadvantage, Coifman, Meyer, and Wickerhauser 
constructed wavelet packets from an MRA wavelet basis. The idea is to decompose the 
wavelet subspaces W 3 into 2° orthogonal components so that each component is generated 
by only the integer translates of a function. In this manner one can get better frequency 
localization. The notion of wavelet packets can be extended to higher dimensions. 


1.7 Organization of the thesis 

This thesis can be divided into two parts. The first part (which consists of Chapters 2-4) 
is on the constructions of band-limited wavelets of L 2 (R) and the Haxdy space H 2 ( R). 
The second part (Chapter 5) is on multiwavelet packets and frame packets of L 2 (R d ) with 
arbitrary dilations. 

In Chapter 2 we construct a set S n for each n > 1 and then we characterize all 
wavelets with Fourier transform supported in S n . The case n — 1 corresponds to the 
Shannon wavelet whereas S 2 is the set associated with the Lemarie- Meyer wavelets. 
Moreover, we show that, if n > 3, then none of the wavelets associated with S n is an 
MRA-wavelet. The case n = 3 also includes the Journe wavelet and the Lemarie wavelet 
presented in §1.3. We have also computed the dimension functions D n associated with 
these wavelets. The functions D n turn out to be even functions and they assume all 
values from 0 to the maximum value attained by D n . In Chapter 3 we construct, with the 
help of the wavelets characterized in Chapter 2, wavelets satisfying interesting properties. 
One of these families answers in positive a question posed in [Gar] about the existence 
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of wavelets with certain properties. We also show that each of the equivalence classes 
M n , n > 0 of wavelets defined in [Web] is non-empty. In addition, a family of wavelets 
in Mo is constructed. Chapter 4 deals with the construction of if 2 - wavelets. We prove a 
result on the structure of interval wavelet sets of H 2 ( R) and characterize all wavelet sets 
of H 2 ( R) which are unions of exactly three intervals. Finally, in Chapter 5 we study the 
multiwavelet packets and frame packets of L 2 (R d ) where we have taken arbitrary dilation 
matrices and allowed the resolution spaces to be spanned by several scaling functions. 



Chapter 2 


A class of band- limited wavelets of 

L 2 (R) 


In this chapter we will construct a class of band-limited wavelets of L 2 (R). For each 
integer n > 1 we construct a set S n and characterize all wavelets whose Fourier transform 
is supported in S n . Then we show that these wavelets are not associated with any MRA. 

2.1 Introduction 

A function is said to be band-limited if its Fourier transform has compact support. The 
simplest example of a band-limited wavelet is the Shannon wavelet whose Fourier trans- 
form is the characteristic function of the set S = [— 2-k, — 7r] U [7r, 2n]. Lemarie and Meyer 
([Meyl], [LM1]) constructed a wavelet whose Fourier transform is infinitely differentiable. 
This wavelet satisfies V>(£) = e^/ 2 &(£), where b is an even, non-negative “bell-shaped” 
function and its support is equal to the set [ — f^r, — |] U [§7r, § 7 r]. In [BSW] the authors 
characterized all wavelets ip such that \ip\ is an even continuous function whose support is 
equal to the set [— 2%— d, — 7r+e]u[7r— e, 27r+e'], where e, e' > 0 and e+e' < tt, so that this 
class includes the Lemarie-Meyer wavelet as a particular case. Later, in [HWW1] (see also 
Theorem 4.1, Chapter 3 of [HW]), the restriction on ip (i.e., \ip\ is even and continuous) 


15 
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was removed and all wavelets with Fourier transform supported in [— § 7 r, — |] U [§ 7 r, |tt] 
were characterized 1 . 

A special case of band-limited wavelets, the MSF wavelets, were considered in [HKLS] 
in which these were called the unimodular wavelets. Recall that for an MSF wavelet 
ip, there exists a measurable set K of R such that |V>| = xk- Denote K + = K n R + 
and K~ = ADR - . The authors of [HKLS] characterized all MSF wavelets for which 
the associated wavelet set K is symmetric with respect to 0 and K + is a union of two 
disjoint intervals. They proved that these are precisely those wavelets for which K + = 
[|^y 7T, 7r] U [2 n 7r, yTif 7r ] for some n > 1. This includes the Journe wavelet. 

This chapter is organized as follows. In §2.2 we will construct a set S n for each n > 1, 
which is symmetric with respect to the origin and is a union of four intervals. In §2.3 
all wavelets ip with Fourier transform supported in S n will be characterized. Next, in 
§2.4 we show that none of these wavelets is associated with an MRA. For a wavelet to be 
associated with an MRA, it is necessary and sufficient that the corresponding dimension 
function is equal to 1 a.e. In §2.5 we will compute the dimension function explicitly and 
show that it is an even function and takes all integral values from 0 to the maximum 
value it attains. This will provide an alternative proof of the fact that the wavelets are 
not MRA-wavelets. 


2.2 The set S n 

Let n > 1. Put 

on—l 

tin = 2 n — 1 
Cn — 2 n — 1 ^ n ^ Qn 2 n — 1 

Define S n = 5+ U S~, where 5+ = [a n , b n ] U [c n ,d n \ and S~ = ~(S+). 

Remark 2.1. 1. If n = 1, then a x = 7 r, bi = c x = d x = 27t and Si is the support of 

the Fourier transform of the Shannon wavelet. 


b n = 2a„ 


-^-7T 
2 n —l J 

o2n-l 


: In fact, they characterized all wavelets with Fourier transform supported in the set [-fa, 4k- f a], 0 < 
a < 7 r. It was also shown that for such a wavelet ip, it is necessary that ip = 0 on [—fa, 27r — fa]. 
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2. For n — 2, a 2 = |7r, 62 = c 2 = |7r and d 2 = |7r. So S2 is the set associated with the 
Lemarie-Meyer wavelets. 

3. As n increases, the set [a n , 6 n ] moves closer to 0 whereas [c n , d n \ moves farther away. 
Also observe that the measure of the set S n is 2(|j^7r + ^rr 71- ) which approaches 
37r as n tends to infinity. 

All terms in the equations (1.3)— (1.6) that characterize all wavelets of L 2 ( R) involve 
dilations by powers of 2 and translations by integral multiples of 27r. In the following 
lemma, for various subsets of S n , we identify those translates and dilates which are possibly 
in S n . We present it in a tabular form for easy reference. 

Lemma 2.1. For n > 3, let S n be as above, and let e n = Then 


£ € the set 

£ + 2 kn £ S n unless 

2 j £ 0 S n unless 

\cL n , e n ] 

k = 0 

3 = 0 

[ 6 n , b n ] 

k = 0 or -1 

j = 0 or n-1 

[c^, d n ] 

k = 0 or — 2 n ~ 1 

3 = 0 or -(n-1) 

[ — j ”^n] 

k = 0 

3 = 0 

[— 6 nj —e n ] 

k = 0 or 1 

j = 0 or n-1 


k = 0 or 2 n ~ l 

j = 0 or -(n-1) 


Also observe that 

[e n , b n ] - 27r = [~b n , -en], ( 2 . 1 ) 

2 = [Cn>^n]> 


and 


[Cn, d n ] ~ 2 • 2” 1 ‘K = [-d n , — Cn]- 
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2.3 The characterization 

The following theorem characterizes all wavelets ip such that supp ij> is contained in S n . 
The case n = 2 is Theorem 4.1 of Chapter 3 in [HW] in which case the wavelets are MRA 
wavelets. 

Theorem 2.1. Let n > 2, ip e L 2 ( R) , supp |CS„ and 6(0 = |^(£)|. Then ip is a 
wavelet for L 2 (R) if and only if 

(i) 6(0 1 for a.e. £ €• [un> &n] hJ [ e n , o n ] 

(%%) 6 2 (0 + 6 2 (2 n_1 0 = 1 for a.e. f € [e n , b n ] 

(tii) b 2 (0 + b 2 (£ - 2k) = 1 for a.e. f £ [e n , b n \ 

(iv) 6(0 = b ( 2 ^ 1 (e - 2k)) for a.e. { € [e n , b n ] 

(v) ip ( 0 = e* e ^ 6 ( 0 , where 9 satisfies 

0(0 + 0 ( 2 n ~ l ( t ; - 2k)) - 0(( - 2k) - 9(2 n ~ 1 0 = (2m(0 + 1) tt, 

for some m(0 6 Z, for a.e. £ 6 [e n , b n ] fl ( supp b) n ( 5^1 supp b). 

Moreover, for n> 3, none of these wavelets arises from an MRA. 

Remark 2.2. 1 . Observe that it follows from the above theorem that \ip\ is completely 

determined by its values on [e„, £>„]. On this set let b = \xp\ be an arbitrary function 
taking values between 0 and 1. Using the relations following Lemma 2.1, \ip\ can 
be extented to other sets of S n with the help of properties ( i ) - (iv). Properties 
(ii), (in) and (iv) can also be written as 

b 2 (0 + b 2 ( = 1 a.e. on [c*, d n ] 
b 2 (0 + b 2 (t + 2n) = 1 a.e. on [- 6 „, -e n ] 

b + 27r ) = b @ a,e ' 0n _Cn ^ 


( 2 . 2 ) 

(2.3) 

(2.4) 
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2. If ( supp b) H (yt-isupp b ) has an empty interior in [e n , b n ], then 6 can be chosen to 
be any measurable function. In particular, we can take <9(£) = 0. 

To prove Theorem 2.1 we will use Theorem 1.1 which characterizes all wavelets of 
L 2 (R). 

Proof of Theorem 2.1 

First let us assume that ^ is a wavelet for L 2 (R), supp $ C S n and b = \^\. So ip 
satisfies equations (1.3)— (1.6). Consider equality (1.5) : 

y: 6 2 (£ + 2k*) = 1 for a.e. £ G R. 

kez 

We use Lemma 2.1 to pick out the non-zero terms in this sum. If |£| G [a n ,e n ], then 
only k = 0 will contribute to the sum. So we get &(£) = 1 a.e., which is (i) of the 
theorem. For £ G [e n , b n ], we get non-zero contributions from k = 0 and k — - 1 . That is, 
6 2 (£)+6 2 (£— 27r) = 1 for a.e. £ which proves (Hi). Also if £ € [e n , b n ], then 2 n-1 £ G [c n , d n ], 
and we get (non-zero terms now correspond to k = 0 and k = — 2 n_1 ) 

£> 2 (2 n-1 £) + 6 2 (2 n_1 £ - 2V) = 1 for a.e. £ G [e n , b n }. (2.5) 

This can also be written as 

& 2 (£) + Z> 2 (£ - 2 b tt) = 1 for a.e. £ G [cn, d n ]. (2.6) 

Now consider equality (1.3): 

^b 2 ( 2 J ‘£) = 1 for a.e. £ G R. 

For £ G [e n ,b n ], the only possible non-zero terms correspond to j = 0 and j = n — 1. 
So we conclude Z> 2 (£) + 6 2 (2 n_1 £) = 1. This proves (n). Combining (ti) and (2.5) we get 
condition ( iv ). It remains to prove condition (v). For this purpose we consider equality 
(1.6) with j = n — 1 : 


VK£ + 2krc)ip (2 n_1 (£ + 2k*)) = 0 for a.e. £ G R. 
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If £ € [e n , b n ], the non-zero terms in the sum come from k = 0 and —1. Hence, 

^(0^(2 n-1 £) + ^ (£ - 2 tt)^ - 2tt» = 0 for a.e. £ G [e n , 6 n ]. (2.7) 

Using ( lii ), (2.5) and (2.7), we see that for almost every £ G [e n ,6 n ], the vectors 

2ir)) and (^(2"-‘5),0(2”- 1 ({ - 2*))) 

are orthonormal in C 2 . If we let ^(£) = e l6( ^&(£), then it follows that, for some real-valued 
measurable function a, 


e ia ® (e te ®b{Z),e' e tt- 2ir) b(Z - 2 tt )) 

= (-e-“ , ( 2 "- 1 «- 2 '»H2”- I K-27r)),e-' s ‘ 2 "' ,£ >6(2”- 1 f)) , 


for a.e. £ G [e n ,6„]. But from (n), (m) and (zt;) , we know that 

b( 2 n ~ 1 £) = 6(£ - 2 tt) 

m = &(2 n - i (^-27r)) 


for a.e £ G [e n , . 


So (2.8) can be written as 

+ e _l9 ( 2n_1(?_27r) )l 6(0 = 0 ] 

L T J \ for a.e. £ G [e n ,6 n ]. 

e *a(0 c W(C-2ir) _ e -»fl(2»- 1 f) 6(2 n “ 1 £) = 0 I 

This shows that 

e -ta(£) _ e i[^(?)+^(2 n_1 (?-27r))+’r] a e . on je n , 6 n ] fl SUpp b. 

e -ia(0 _ e 2[e(f-27r)+e(2 n - 1 0] a.e. on [e n , 6 n ] n (i^supp b). 

Hence, for almost every £ G [e n , 6 n ] D supp b fl (j^r supp b), we have 

0(0 + e (2 n-1 (£ - 2 tt)) - 0(£ - 2tt) - 0(2 n " 1 £) = (2m(£) + 1) tt, 


(2.8) 


(2.9) 


for some integer- valued measurable function m. This proves (u). 

We now prove the converse. Suppose ip G L 2 (R), supp ip C S n and the function 
6(£) = |$(£)| satisfies conditions (*) - (v) of the theorem. By Theorem 1.1 it is sufficient 
to show that ||^|| 2 = 1 and ip satisfies (1.3) and (1.4). We have, 

2xWl|| = IW-llt = [ »m 2 dt; = 
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a c n P 6 n P"~0>n P^n pb n pdn \ 

+ / + / +/+/+/ )^(Ode 

■a n J —b n '/-Cn */a n ** 671 J Cn / 

/ -C n /»d n P~~Gn pb n 

+ + / + / +2(e n - a n ). 

■dn J —bn J tn 


By changing variables f -> 2 n_1 (f - 2 tt),£ -+ 2 n ~ 1 f and £ - 2 tt in the first, second 

and third integrals respectively, we get 




1 / [6 2 (2”“ 1 (f - 2x)) + i. 2 (2"- 1 f)] if 

Je n 


P0 n 

+ / [6 2 « - 2tt) + i> 2 «)] dl, + 2(e„ - a.) 

Je n 


2 (b n e n ) + (6 n e n ) -t- 2(e n o n ) — 27 t, 


where we have used (2.5) and property (in) of the theorem. Hence, ||^|| 2 = 1. 

We will now show that ip satisfies (1.3). Let 

p(?) = £ | fta'f) | 2 . 

jez 

Suppose £ > 0. Observe that 

R + = |J 2 l [a n ,b n ] 

le z 

= i_J ^ ([o n , e n ] U \e n , 5 n ]) . 

le z 

So there is an Z € Z such that £ € 2*[a„,e n ] U 2*[e n ,Z> n ]. If £ 6 2*[a n ,e n ] then since 
2“*£ e [ a n ,e n ], by Lemma 2.1 we have, 2 J (2 -J £) fi S n if j ± 0. That is, 2 J £ S n if 
j A -1- Hence, p(£) = 1, by (i). Similarly by using (ii), we can prove that p(£) = 1 if 
£ € 2 i [e n , b n ], A similar decomposition for £ < 0 proves that p(£) = 1 for a.e. (el 
Finally, we have to show that ip satisfies (1.4). For m € 2Z + 1, let us denote the 
function on the left hand side of (1.4) by i m (£). Then 

t m (() = 53lf(2 , {)«(»(f + 2mx)) 

J>0 

— ^ ip (2P (£ + 2m7r - 2m7r)) ip (2? (£ + 2m7r)) 

j>0 

= tlj^+2rrvn). 
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Therefore, we have only to show that t m (£) = 0 for a.e. £, if me 2Z_ + 1 (negative odd 
integers). Suppose m € 2Z_ + 1 and m ^ —1. Let £ € K and suppose that 2 J £ € S n . 
Since j > 0, n > 2 and m/-l and odd, therefore 2 J m ± 0, ±l,±2 n_1 . So by Lemma 
2.1, we have 2 J £ + 2 • 2 3 vmi (f S n . This implies that each term of t m (£) is zero which proves 
that t m is zero. It now remains to show that i_ x (£) = 0 for a.e £. We have 

mo = E^(2 J oS 2j e-2-2^). 

J>0 

Let £ € R and suppose that 2^£ € 5„. Then again by Lemma 2.1, 2>£ - 2 • 2% £ S„ if 
-2 J A 0) ±1) ±2 n_1 . So the only possible f s to contribute a non-zero term are j = 0 and 
j = n — 1. Thus, 

t-i (£) = $(£)$(£ ~ 2w) + ^(2 n - 1 ^)V>(2 n - 1 e - 2 • 2"-%). (2.10) 

Now, both £ and £ — 2-7T belong to S n only if £ E [e n , 6 n ]. Hence, the first term of (2.10) 
is zero unless £ € [e n , b n \. Similarly, both 2 n_1 £ and 2 n_1 £ - 2 • 2 n ~ 1 -jr belong to S n only 
when 2 n-1 £ 6 [c n , d n ], which is equivalent to saying that £ E [e n , 6 n ]. That is, the second 
term of (2.10) is also zero unless £ € [e n , £>„]. Thus we get, £_i(£) = 0 if £ £ [e„, 6 n ]. Now 
on [e n ,6 n ], we have, by (2.9) 

t - i {0 = K £) K 2 n " 1 £)^ (f) " s({ ' 2,)1 + H 2 fl ' 1 £) K ^ !l<(2 ’‘ 1?)+ ^ 2 " 

If £ 0 [e n , 6 n ] n supp 6n (jT^rSupp 6), then either &(£) = 0 or 6(2 n_1 £) = 0. So t_ x (£) = 0. 
And if £ € [e n , b n ] n supp b n (^r supp b ), then by (v), i_ x (£) = 0. This completes the 
characterization. □ 

There is another equation which characterizes all wavelets ip such that ip is supported 
in the set S n and the function b = \ip\ is even. For n = 2, the following proposition 
is proved in [HW] (see Proposition 4.7, Chapter 3). We observe that the result can be 
extended to the general case. 

Proposition 2.1. Suppose that ip is a wavelet o/L 2 (M), b = \ip\, and supp b C S n . Then 
b is almost everywhere even if and only if 

& 2 (£) + b 2 ( 2n - £) = 1 for a.e. £ € [e n , 6 n ]. 


( 2 . 11 ) 
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Proof: Suppose that b is an even function and f G [e n , b n ] Since -£ G — e„] , we 

have 

1 = 5 2 (-0 + & 2 R + 27r), 5y(23 ) 

= 6 2 (0 + 6 2 (27r - £), since b is even, 

which is (2 11) 

Conversely, suppose that (2 11) holds Since by (i) of Theorem 2 1, b is even on [a n , e„], 
it is enough to show that b is even on the sets [e n , b n ] and [c„, d n ] 

(a) Let f € [e n ,b n ] Therefore, — £ G [—b n , —e n ] Then 

b 2 (- 0 + 6 2 (-{ + 2ir) = 1, by (2 3) 

This fact, together with (2 11), gives us 

m = K-fl ( 2 12) 

(b) Now let £ e [c„, d„] Therefore, £ [en, f„] FYom (2 12) we get 

(213) 

Since - € [-& n , -e„], using (2 3) we obtain 

^(-2^re)+& 2 (-2^ + 2*)=l (2 14) 

Now, since -£ G [-d n , — c„], we get (usmg (2 4)) 

»(-^I« + 2rr)=K-0 (2 15) 

Substituting (2 13) and (2 15) in (2 14), we get 

62 (2^Te) +6 2 (-« = l. (2 16) 

Comparing (2 16) and (2 2) we get &(£) = 6(-£) This proves the proposition □ 
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2.4 The wavelets associated with S n , n > 3 are non- 
MRA 


We will now show that, if n > 3, none of the wavelets characterized in Theorem 2 1 is 
associated with an MRA 

Let ip be any wavelet such that ip is supported in S n and let it be associated with an 
MRA Let cp and mo be the corresponding scaling function and low pass filter respectively 
Then we have ( see equations (1 11) and (1 7)) 


i^(«i 2 =£h&(n)i ! 

j>i 


and 


$(2£) = mo(Ofe) 

Using (2 17) we can easily find that 

t 

1 if |£| < 0„ 

\k0\ = •{ b(2 n ~ l ~ 1 0 if |£| e 2 z [e n , b n ], 0<l<n-2 
0 otherwise 


That is, 


if iei < f 


\tm = < 


b( 2 n ~ l £) if |£| € 2 z_1 [e n ,& n ], 0 < / < n - 2 
0 otherwise 

Case 1: &(£) ^ 1 on [ e n ,b n ], i.e , &(£) ^ 0 on [c„,dn3 
If |£| < On> then from (2 18) and (2 19) we have 


(2 17) 


(2 18) 


(2 19) 


(2 20 ) 


|$(2£)l = l m o(OI lfe)l = l™o(£)l 


Therefore, by (2 20), we obtain 





A class of band-limited wavelets of L 2 ( R) 


25 


Since m 0 is 27r-penodic, we have 


K(OI = 1 


on 


fl'n Ojti 

r 


+ 2 2 n— 3 7r 


Note that [-f , f] + 2 2"- 3 7r = [f, f] = 2 n " 2 [e n , b n \ Now, from (2 19), on 2 n - 2 [e re , b n ] 

we have |0(O| = 6(2^ So |0(2O| = K(0I |0(OI = 6(20 But by (2 20), |0(2O| = 0 

on 2 n-2 [e n , 6 n ] Therefore, 6(20 = 0 on 2 n ~ 2 [e n ,6 n ] That is, 6(0 = 0 on 2 n ~ 1 [e n ,6 n ] = 
* 

[cn, d n ], which is a contradiction 
Case 2: 6(0 = 1 on [e n ,6 n ] 

As above, we have 


Therefore, 


1 if |£|<a» 


m\ = { 


or if f € 2 l [-6 n , — e n ], 0 < l < n - 2 


^ 0 otherwise 


(2 21 ) 


1 if |0 < * 

10(20 1 == ^ or if £ G 2 i-1 [-6„, -e n ], 0</<n — 2 
0 otherwise 


(2 22 ) 


Now |0(2OI = ho(OI |0(OI = ho(OI if 10 < Qn, by (2 21) Therefore, by (2 22), 
|m 0 (OI = 1 if |0 < Using the 27r-penodicity of m 0 , we get 


ho (01 = 1 on [-y>y 


2 2 n " 3 7T 


Note that 


On 2 n-2 [— 6 n , — e n ], we have 


10(201 = K(0I ltel = K(OI = 1 

But by (2 22), |0(2O| = 0 on [— ^f,— ^], which is again a contradiction Therefore, the 
wavelets of Theorem 2 1 are non-MRA wavelets, if n > 3 
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2.5 The dimension functions of wavelets associated 
with S n 


Recall that for a wavelet ip of L 2 ( R) to be associated with an MRA, it is necessary and 
sufficient that — 1 a e , where 


i V«) = £EhH 2, tt+ 2 *’0) P 

3> 1 k€ Z 

In this section we will compute the dimension functions for the wavelets characterized 
in Theorem 2 1 This will provide an alternative proof of the fact that these wavelets are 
not associated with any MRA These dimension functions are symmetric with respect to 
the origin and attain all integral values starting from 0 to the maximum value Also note 
that if 77 . > 3, then all wavelets whose Fourier transform is supported in S n have the same 
dimension function D n We will prove the following proposition 


Proposition 2.2. Fix n > 3 Let ip be a wavelet such that supp ip C S n and let D n be 
the dimension function associated with ip Then 


n — 1 
r — 1 


Dn( 0 = < 


0 

1 


a e. if |f | G [0, ^tt] 

a e if |f | € (2 < r < n - 1) 

a e if | G Jrrf ?r] = [a n , e n ] 

a e if |f | G [pzf7r,7r] = [e n , 7r] 


Proof: For l G Z, we define 


2*1+1 — i 2 n+i_1 2 l 

Pi = ^ — 7 * = and q t = — ^ — 7 — tt = 2 _1 e n 


2 n - 1 


2 n - 1 


Note that p 0 = a n , pi - 2 p 0 = b ny p n = d n , q x = e n and q n = Cn Also observe that 
Pi < q i+ 1 < p i+ 1 for all Z > 0 As D n is 27r-penodic, it is enough to compute its values for 
f G [ — 7 r, 7 r] Again as [-7r,7r] C [~PuP\], we will compute D n for the interval [-pi,Pi] 


(0,pi] = |J2 l \p 0 ,pi] = |J[p-i, P- 2 + 1 ] 

i >0 />0 


Note that 
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Case 1: £ 6 \po,Pi] 

If 2 m < k < 2 m+1 — 1 (0 < m < n - 3) then £ + 2kn G [p n+ %, g m+ 3 ] So 2 J (£ + 2/br ) G 
bj+m+ 2 , 9j+m+3] Ifj+m + 3<n then 2 J (£ + 2fc7r) G [pi> (7n] = [4, Cn ] which is not m 
supp ip, and ifj+m + 3>n+l then 2 3 (£ + 2kir) > d n , which shows that 2 J (£ + 2kn) £ 
supp 'ip Now if k > 2 n ~ 2 , then £ + 2kir > d n Thus, we have proved that if k > 1, 
then 2 J (£ + 2kix) £ supp ip for all j > 1 Similarly, we can show that if fc < —2, then 
2 J (£ 4- 2kir) supp -ip for all j > 1 Therefore, we have only to consider k = 0, -1 Note 
that bo, Pi] = bo. ft] U [ft.pi] 

(a) If £ € bo> ft], then 2 J £ € fa, ft +1 ] Now j > n =► 2 J £ > p n = d n and j < n - 1 =4 
J + 1 < n =>- for all j >1, 2 J £ G fa, q n ] C [pi, ft] = [6 n , ft ] which is not in supp ip Also, 
£ — 27r G [—ft, — Pi] =>• 2 J (£ — 27r) G [— g J+ 2 , — p 5 +i] A similar argument as m the case of 
2 J £ shows that 2 J (£ - 2n) g supp ip for all j > 1 So !>„(£) = 0 a e on bo. Pi] 

(b) If £ G [ft.pi], then 2 J £ G [ft + i,p J+1 ] So 2 J £ # supp ip if ; ^ n - 1 Similarly, 
2 J (£ - 27 t) G [— Pj+i, — ?j+i] So 2 j (£ - 2i r) ^ supp ^ if y ^ n — 1 Therefore, D n (£) = 
|'0(2 n-1 £)| 2 + | ip (2”- 1 (£ - 2k)) I 2 = 1, by (2 6), as 2 n " 1 £ G [q n ,p n ] = [ft, 4] 

We now proceed to compute D n on other sets Observe that 

Dn(t) = J2 Y^\$( 23 tt +2k *))\ 2+ Yl$( 2J &\ 2 

k=—oo j>l ]>1 

+ EEI^ 2) k+ 2 ^))I ! 

A:=l J>1 

= D n (0 +^°(0 + ^n ( 0 . sa y 

Let £ G = b-i+i.P-i+ 2 ], l > 2 

As in the above cases, we can show that if k < — 2 n ~ 2 , or if -2 m+1 -1 < k < — 2 m — 1 
(0 < m < n - 3), then 2 J (£ + 2 /ot) £ supp ip for all j > 1 Also, if A = -2 m (0 <m < 
n - 3), then 2* (£ + 2fc7r) £ supp ip for all y > 1 unless j = n-m- 2 Hence, m the sum 
D~(£), we have only to consider j = n — m — 2 (0 <m < n — 3) and k — —2 m Thus, 

B ;(«= E |’?( 2 "' m ~ 2 «- 2 2m,r ))f 
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In a similar manner, we can show that 

71-3 ^ 

oj«) = E|^( 2 ”"”“ 2 tt+ 2 *"*)) 

771=0 

Case 2: £ G [|^;7r, 5zy-7r], Z > n 

By a straightforward calculation, one can show that if l > n, then 

2”— m- 2 (£ + 2 2 m 7r) € [c„, d n ], for 0 < m < n — 3 

Therefore, 

^(2 n " m - 2 (^ + 2 2 m 7r)) + ^(2 n_m-2 (^ — 2 2 m 7r)) 

= 2 n-m-2(^ + 2 2 m 7r)) +^(2 n - m " 2 (e + 2 2 m 7r) - 2 n - 1 7r) 

= 1, by (2 6) 

Hence, we get D+ (£) + £>"(£)= n - 2 Now, [p_z+i,p_z +2 ] = [p-z+i, q-1+2] U [q- l+2 ,P-i+ 2] 

(a) If £ € [p-z+i, ?-z +2 ], then 2 J £ £ supp if j ^ - 1 But 2' _1 £ e [p 0) ?i] = [oi, e n ] 
Therefore, £>£(£) = |^(2 i " 1 0l 2 = 1 

(b) If £ e [q-i + 2 , P- 1 + 2 ], then 2 J £ & supp ^ if .7 7^ Z — l,n H- / — 2 Now 2’ l+!-2 £ 6 
[?n,Pn] = [Cn,dn] By (2 2) we get £>£(£) = |^(2 ?l+ *- 2 £)| 2 + |V’(2 f_1 0l 2 = 1 Thus, 
in either case D n {f) = n — 2 + 1 = n — 1 We have proved that D n (£) = n — 1, if 
ee [g^r,^], Z > n 

Case 3: £ e [f^Tr, ?^r 7r ] » 1 = 2 

Again as in Case 2, it can be shown that 

2 n— m— 2(£ _ 2 2 m 7r) and 2 n-m-2 (£ + 2 2 m 7r) 

are not m supp for all m such that 0 < m < n — 3 Therefore, 

£>;;(« = = o 

Now if £ € [|^7r, 2 2»li 1 ?r ] ^ en 23 f £ SU PP $ if J J 4 1, and if £ S [^nrr 7 !, f^ir] then 
2+ £ supp if j 1, n In either case, 


A.(£) = Bj(f) = 1 
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Case 4: 3 < l < n — 1 (This case is required if n > 4 ) 

We can write 2 n ~ m ~ 2 £ G tt, f^7r] , where p = n — m — l — l Using the condition 

on support of ip, we can show that if £ € , then f ± 2 2 n ~ 2 7T e supp ip 

only when p < -1 So m order that ^(2 n_m_2 ^ ± 2 2 n_2 7r) ^ 0, we must have p = 
n — m — l — l < — 1 , i e , m > n — l We also have, 0 < m < n — 3 Therefore, 

1 < n — m — 2 < l — 2 Thus 

A,© = ^|-A(2 J f)| 2 

J>1 
n— 3 

+ E{l^( 2n “ m ' 2 (£ + 2 2 m 7r))| 2 + |^(2 7l " m - 2 (e-2 2 m 7r))| 2 } 
m=0 

= El^>f + E{l^ + 2 2”- 2 7r)| 2 +|^(2'5-2 2 n-2 7r)| 2 } 

J>1 3=1 

Now, 2 s £ + 2 2 n ~h r G [cn, d n ], 1 < s < Z — 2 By (2 6), the second sum on the right hand 
side of the last equality above is equal to l — 2 Also as in the previous cases, the first 
sum can be shown to be equal to 1 Hence, D n (£) = 1 — 1. □ 

2.6 Concluding remarks 

In [HKLS] MSF wavelets were considered where they were called ummodular wavelets 
Let ip be an MSF wavelet and K = supp ip Define K + — K n R + and K~ = K fl R - 
One of the results proved m [HKLS] is about wavelets ip such that K~ = —K + and K + 
consists of two disjoint intervals They proved the following 

(A) Let ip(£) = p{Oxk{£)i where K = K + U K~,K~ = -K + and K + = [a„, ix) U 
[2 n-1 7r, d n \ with \p(()\ = 1 Then ^ is a wavelet for L 2 (R) Moreover, each MSF wavelet 

i 

for which K~ = -K + , and K + is a union of two disjoint intervals is of this form 

If we define b = X[e n ,A on [e„, b n \, and extend it to S n by using (i)-(iv) of Theorem 

2 1, then we get the wavelets described in (A) In particular, for n = 3, the corresponding 
wavelet is the Journe wavelet 

The paper [HKLS] also contains the following result 
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(B) Let V>(£) = h{Z)xk{£) where 


K~ = 



2p 

2 n 


+ 1\ /. 2p+l\ ' 


and 


K + = 

'2{p + 1) 2(2p + l) ‘ 

u 

’2 n (2p + 1) 2 n+1 (p+l) 

7T, 7 r 


2 n - 1 ’ 2 n - 1 


2"-l ’ 2 n — 1 


for n > 3, 1 < p < 2 n ~ 1 -2 and |p(£)| = 1 Then if> is a wavelet for L 2 ( R) Moreover, each 
MSF wavelet for which K~ is an interval and K + is the union of two disjoint intervals is 
of this form 

If we define 6(£) = 0 on [e n ,f> n ] and extend it to S n by using (i)-(iv) of Theorem 2 1, 
then we get 

^ X( — &n, - 'fln]CJ[an)Cn]U[cn,dn] 


which corresponds to the case p = 2” 2 - 1 in (B) In particular, for n = 3, we get the 
Leman6 wavelet 



Chapter 3 


Some more band-limited wavelets of 

L 2 (M) 


In this chapter we will construct wavelets with interesting properties These wavelets will 
be generated with the help of the wavelets characterized m Chapter 2 We will construct 
two families of band-limited wavelets with Fourier transform discontinuous at the origin 
One of these families have Fourier transforms which are even functions We will also show 
that the equivalence classes of wavelets defined in [Web] are non-empty 

3.1 Wavelets with discontinuous Fourier transform 
at the origin 

Every wavelet ip of L 2 (R) has to satisfy (1 3) Therefore, ^(0) must be equal to 0 if \ip\ 
is continuous at 0 In fact, if ip is a band-limited wavelet and \ip\ is continuous at 0, 
then a stronger result holds |$| = 0 a e m an open neighbourhood of the origin (see 
[HW, Theorem 2 7, Chapter 3]) In [HWW1], it was shown that if ip is a wavelet with 
supp ip c [-47T + fa, fa], 0 < a < 7r, then \ip\ = 0 on the interval [-27T -I- fa, fa] It 
is interesting to observe the behaviour of $ at the origin if we enlarge the support of 
Some related results are the following 


Ol 
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(1) Madych [Mad] constructed an example of a band-limited scaling function so that 
the corresponding wavelet ip (which is also band-limited) has the property that ip does 
not vanish in any neighbourhood of the origin 

(2) Garngos [Gar], m his Ph D thesis, constructed a wavelet ip with supp ip c [ — 47 r, 7 t] 
and ip not vanishing in any neighbourhood of the origin 

(3) In [BGRW], the authors constructed for each e > 0, a wavelet ip € such that 
supp ip € C [— |7r, |7r + e], but ip £ does not vanish in any neighbourhood of the origin 

Recall that (see Definition 1.2, Chapter 1) a wavelet ip is called a minimally supported 
frequency (MSF) wavelet if \ip\ = xk for some measurable subset K of R. Such a set K 
is called a wavelet set, which necessarily has measure 27r 

In §3 1 1 we will construct a family of wavelets with Fourier transform not vanishing 
in any neighbourhood of the origin The wavelets we are going to construct (as well 
as the wavelets mentioned m (l)-(3) above) are MSF wavelets so that it is sufficient to 
construct the associated wavelet sets For n > 2, we will start with a wavelet set W (which 
depends on n) and construct another set W n Then we will show that W n is both dilation 
and translation equivalent to the wavelet set W, which will prove, by Corollary 1 2 of 
Chapter 1, that W n is also a wavelet set 


3.1.1 A family of band-limited wavelets with Fourier transform 
discontinuous at the origin 


Let n > 2 Put 


Let 


(In — 


HT 

2 n — 1 ' ’ 


b n — 2d n 2 n — 

0. = ^-^, 4 = 2X = &1 


_ 2 n -2 
2 n -l 


7 T 


L = [ -d n , -Cn], M = [~e n , -a n ] and R = [a n , b n ] 


In Theorem 2 1, by taking &(£) = 1 on [e n , b n \ and extending to the whole of S n , we get 
a wavelet whose Fourier transform is the characteristic function of the set L U M U R 
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Therefore, W - LUMUiJisa wavelet set Let e be a positive number such that e < 
Define the following sets 


A 


a n . € a n 

I.T + F’T + e J ’ 


P 2 — [o„ + 2e,Z> n ], 

P 3 = [d n , + 2e] 

To make P 2 a non-empty set we need 2e < a n Let 


X 0 = P 1 -2 2 n ~h r, F 0 = iX 0 , 

Xi = — 2 • 2 n-2 7r, 1^ = 2 ^^, Z > 1 

Now define the set 

00 00 

W n = (L\U^) U (U^) UMU ( p i up 2 u p 3) 

1=0 1=0 

Theorem 3.1. The set W n , defined m (3 1), is a wavelet set forn>2 


( 31 ) 


Proof: We will prove the theorem by showing that W n is translation and dilation equiv- 
alent to the wavelet set W First of all, we show by induction that Xi C L for all Z > 0 
Note that Pi C [-o n , a n ] So 

X 0 = Pi - 2 2 n— 2 7t C [— a„, a n ] - 2 n ~ l ir = [~d n , -c n ] = L 


Now, assume that X m c L Then, 

Ym = 2m+n^ m C [" 


d n 


2m+n 5 2 m + n J 


rn+ n 1 C [ a n> a n] 


X m+ i = Y m - 2 n 1 'K C L 


The interval Xo lies in the interval [— 2 n ~ 1 7r, — c n ] and {Xi l > 1} lie m [—d n ,—2 n ~ 1 ] 
Further, Xi+i lies to the right of Xu l > 1 The intervals Yi, l > 0 he in ^ [— a n , ] so 
that Yi+i lies to the right of Yi for all Z > 0 Also observe that {Xi,Yi l > 0} is a disjoint 
collection 

Dilation equivalence 

2P 1 UP 2 U^P 3 = [a n + 2 ^r,a n + 2e] U[a n + 2e,6 n ]U [a n ,a n + ^r_ 
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— [ci n , b n ] — R, 

OO OO 00 oo 

(i \ U X,) U (U 2" + ‘V',) = (l \ U X,) u ftj x,) = L 
1=0 1=0 1=0 1=0 

(since Xi C L , for all l > 0) The set M appears in both the partitions of W and W n 
Translation equivalence 

P 2 U (P 3 - 2 2 n-2 7r) = [o n + 2e, b n ] U [a n , o„ + 2e] = R, 

00 OO 

(L\y^)u(|J(^-2 2' 1 - 1 tt)) U (Pi - 2 2”- 2 7r) 

1=0 1=0 y 

00 OO 

= (£\U*l)u(U*i)uJf 0 

1=0 V l = 1 

= L 

Again, M appears in both the partitions of W and W n As W n is translation and dilation 
equivalent to the wavelet set W, we have proved that W n is a wavelet set □ 

Let 'ip n be the characteristic function of W n Then ip n is a band-limited wavelet Since 
Yi c ^[-a n , for all l > 0, $ n does not vanish in any neighbourhood of 0 In 
particular, ip n is discontinuous at the origin 

3.1.2 Band-limited wavelets with Fourier transform even and 
discontinuous at the origin 

In [Gar] the question of the existence of a wavelet with the following three properties was 
asked 

(i) ijj is band-limited, ie., $ is compactly supported, 

(ii) ijj is even, and 

(iii) ip does not vanish in any neighbourhood of the origin 



Some more band-limited wavelets of L 2 ( R) 


35 


Examples of wavelets satisfying any two of the above three properties can be constructed 
The Shannon wavelet ip (where ip = satisfies (1) and (n) but not (in) The 

wavelets constructed in §3 1 1, as well as the wavelets of [Mad], [Gar] and [BGRW] referred 
in §3 1 1, satisfy (i) and (m) but not (n) N Arcozzi [Arc] constructed a wavelet which 
satisfies (n) and (m) but not (i) In this section we will construct a family of wavelets 
satisfying all the three properties listed above These wavelets are again MSF wavelets 
and so we construct the associated wavelet sets To get these wavelets we will suitably 
modify the method presented in §3 1 1 

For n > 2, let a n , b n , d n and e n be as in §3 1 1 
Let 

L\ = [-7T, —On], L 2 = [~d n , — 2 n_1 7r], 
and 

Ri = [<*»> tt], #2 = [2 n_1 7r, d n ] 

On [e n , b n ], define b = X[e n ,*] Then, by extending b to the whole of S n (see Theorem 2 1), 
we see that the characteristic function of the set L\ UL 2 URi UR 2 is the Fourier transform 
of a wavelet Therefore, it is a wavelet set For e > 0 such that z < \ (pzf 7r) = ^, 
construct the following sets- 

Si 

S 2 

S * 
and 

7, — — Uj 1UX li — J-, o 

To make 5 2 a non-empty set, we need to take o n + 2e < 7r which is equivalent to e < ^ L 
Let 

Eq ~ Si + 2 2 n— 2 7T , Fo = 

E x = FU + 2 2"- 2 tt, F l = ^ m E l , l> 1, 

Gi = —Ex and Hi = —F t , for l > 0 


■ [ 


a JL + ± + e 

2 2 ”’ 2 . 


= [a n + 2e,7r], 
= [d n , d n T 2e] 


_ 1 O O 
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Define 

00 oo 

Kn = [R 2 \ U^i) u (|J F 0 u (& U S 2 U S z ) 

V 1=0 1=0 
OO OO 

u [L 2 \ U Gi) U (U Hi) U (Ti U T 2 U T 3 ) 

i=o {=0 

Theorem 3.2. For each n> 2, the set K n , defined m (3 2), *5 a wavelet set 


(3 2) 


Proof: We will prove that K n is translation and dilation equivalent to the wavelet set 
Li U L 2 U R\ U i ?2 As in §3 1 1, we can show by induction that E\ C F 2 for all Z > 0 
By symmetry of the set K n ,Gi C L 2 for alH > 0 The intervals Ei, l > 0 lie inside the 
interval [2 n_1 7r, rf n ] and F/ + i lies to the left of Ei for all Z > 0 Similarly, the intervals Fi, 
l > 0 lie in ;jr so that F/ + i lies to the left of F/ for / > 0 Similar statements are 

true for the intervals Gi and Fi, l > 0 
Dilation equivalence 

2 Sx US 2 U 

and 

00 

(ft \QjSj) U (U 2 n+,+1 F,) = (ft \ U E ‘) u (U E ‘) = 

1=0 1=0 1=0 1=0 

(since Ei C R 2 , for all l > 0) Similarly for L x and L 2 
Translation equivalence 


[ a n + 2 n^T’ a n + U [a n + 2e, 7r] U |a n , 0 n + 2 ^rr] 
[a„, 7 r] = Fi 


OO 00 00 ^ 


S 2 U (^3 — 2 2 n 2 7 r) — tz n + 2 €, 7 r] U [a nj a- n + 2e] — Ri 


and 

o° °o . 

(f 2 \QEi)u(|J(F ! + 2 2 n - 2 7r))u(S' 1 + 2 2 n - 2 7r) 
v 1=0 1=0 

00 00 

= (f 2 \ 0F() U (U Fi) U Fo = F 2 
1=0 1=1 

Similarly for Li and L 2 Thus, we have shown that K n is both dilation and translation 
equivalent to the wavelet set L\ U L 2 U R\ U i? 2 , and hence is a wavelet set □ 
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Let w n be the characteristic function of the set K n So w n is a band-limited wavelet 
such that w n is even Since F; c ^r[|, for all l > 0, w n does not vanish in any 
neighbourhood of 0 In particular, it is discontinuous at 0 

3.2 Dimension functions of band- limited wavelets 

Let us recall that for every wavelet ip, there is an associated function D^, called the 
dimension function 

a*(o=EE i^( 2 j «+ 2 ^)) i 2 

3>i fcez 

It is clear that D $ is 27r-penodic and is m L 2 (T) First of all, let us observe that the 
dimension function is bounded for every band-limited wavelet 

Proposition 3.1. If ip is a wavelet and supp ip C [— 2mr, 2wr] for some n 6 N, then 
D^<n 

Proof: Let 

m = D^oi 2 

j>i 

By the condition on supp ip, we have supp F C [-n7r, nir] Also, since |^(2 J £) | 2 = 1 
for every wavelet ip (see (13)), we observe that F < 1 Therefore, F < X[-mr,mr] This 
implies that 

■Dy >(£) = } ^ F(£ + 2H) < ^ ] X[-n 7 r,Tnr](^ -h 2k-K) = n, 

kez fcez 

which proves the proposition D 

The statement of the above proposition is not sharp It was proved in [HW] (see also 
[HWW1]) that if supp ip C [~§7r, |tt] , then ip is an MRA wavelet, ie,^ = lae Now 
consider the MSF wavelet constructed in §3 1 1 for the case n = 2 The corresponding 
wavelet set is W 2) as defined m (3 1) If £ € [|7r, |7 t + e], then 

r 8 8 

4(^-27r)e — -7T, --7T + 4e 


“8 g 

and 2£ G -tt + 2e 
o o 
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Both the sets [— §*"} — |vr + 4e] and [§7r, |7r + 2e] are contained in the wavelet set W2 
Therefore, on the set [|7r, |7r 4- e], the dimension function is at least 2 (the fact that it is 
at most 2 will follow from the proposition we prove below) So by increasing the support 
of ip slightly larger than [— |7r, |7 t], we no longer have the dimension function of ip bounded 
by 1 That is, [— §7r, §7r] is the largest symmetric interval to conclude that the dimension 
function is bounded by 1 whenever ip is supported in this interval This motivates to find 
out other symmetric intervals of the form [~t n ,t n ],n > 2 such that supp ip C [~t n ,t n ] 
implies that D $ is bounded by n The following result provides a lower bound for t n 

Proposition 3.2. If ip is a wavelet and supp ip C [— *•] for some n € N, then 

D-,1, <n 

Proof: Since Dy is 27r-penodic, we will prove that if ip satisfies the hypothesis, then 
D^,(£) < n for £ € [ — 7r, 7r] For £ € [ — 7r, 7r], we have 2 3 (2k — 1)tt < 2 J (£ + 2kir) < 
2 3 (2k + l)n 

( I ) 3 = n If k > 2, then 2>(£ + 2A:7r) = 2 n (£ + 2 Ajtt) > 2 n (2k - l)ir > 3 2 n 7r > 
Similarly, if k < -2, then 2 n (£ + 2kn) < ~^ 7 r Hence, for j = n, the only non-zero 
terms contributing to are for k = —1, 0, 1 

( II ) 3 > n + 1 If k > 1, then 2 3 {£ + 2kit) > 2 3 {2k - 1)tt > 2 n+1 (2k - 1)tt > 2 n+1 7r > 
Similarly, if k < —1, then 2 3 (£ + 2kn) < ^p7r Hence, for j > n + 1, only 

contributing k to is k = 0 
Thus, we have 

= EE I ^(2 J K + 2*»))| 2 

]=i fee z 

+ E Il5(2"« + 2M)I 2 + E l^ 25 «l 2 < 33 ) 

A=— 1 J>n+1 

= EEi^ +2 **))i 2 

j=i ke z 

+ { \ip (2«(e - 2 tt)) i 2 + $ (2 n (e + 2*0) 1 2 } + E i 2 ( 3 4 ) 
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If £ € [-§*, §7r], then 2"(C + 2 tt) > 2 n |tt = ^tt Similarly, 2 n (£ - 2 tt) < -^tt 
S o both the terms in the second summand in the right hand side of (3 4) are zero and 
D+ <n-l + l = n, where we have used equations (1 3) and (1 5) Now, if £ € [§ 7 r, 7 r], 
then for all j > n + 1, 2 J £ > 2 n+1 |tt = So the last sum in (3 3) is zero and 

again <n In a similar manner, it can be shown that A/, < n if £ e [- 7 r, — § 7 r] This 
finishes the proof of the proposition □ 

Proposition 3 2 was also proved by Z Rzeszotmk and D Speegle in an unpublished 
article In the same article they also proved that for every n > 2 and 0 < e < S(n), there 
exists a wavelet %j> with Fourier transform supported in the set [-^ 7 r, ^ 7 t + e] such 
that D^> n on a set of positive measure This shows that is the largest 

symmetric interval to conclude that the dimension function Df is bounded by n whenever 
is supported in this interval We are grateful to Prof Guido Weiss for providing this 
information to us 


3.3 Equivalence classes of wavelets 

First of all, we will recall the equivalence relation on the set of all wavelets of L 2 (R) 
defined in [Web] Let ijj be a wavelet of L 2 (R) Define V 3 — sp{ipi,k l < J, k € Z} Then, 
it is easy to verify that the subspaces V v j € Z satisfy the following properties 

1 V 3 C V J+ i for all j e Z 

2 / € Vj if and only if /( 2 ) 6 V^+i for all j G % 

3 Uj gzVj = L 2 (R) and fljgz = {0} 

4 Vo is invariant under translation by integers 

A natural question is Are there other groups of translations under which Vo is invariant 7 
For ol 6 R, let T a denote the unitary operator T a f = /( — a) Consider the group of 
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translations Q n = { T. j » m € Z} and the group — {T a a 6 i} Let C n be the 
collection of all wavelets such that the corresponding space Vo is invariant under the 
group Q n Then we have 


£o D Ci D C 2 D • C n D C n+ 1 O ■ D Coo 

An equivalence relation can now be defined on the collection of wavelets The equivalence 
classes are given by M n = C n \ C n+ 1 , with Moo = £oo Therefore, M n is the class of 
wavelets such that Vo is invariant under the group Q n but not under Q n+ \ 

In [Web] Weber proved the following facts 

(l) Ml oq — C\n>QCn 

(u) Moo is precisely the collection of all MSF wavelets 

For the classes M n ,n > 1 the following characterization is given 

(ui) The equivalence class M n consists of all wavelets ip such that supp ip is not partially 
self similar with respect to any odd multiple of 2*7 t , k = 1,2, ,n but is partially 
self similar with respect to some odd multiple of 2 n+1 7T 

Definition 3.1. [Web] A set E is said to be partially self similar with respect to Oi 6 M if 
there is a set F of positive measure such that both F and F + oi are subsets of E 

s In the same paper, Weber produced examples of wavelets belonging to the first few 
equivalence classes Mm namely for n = 0,1,2 and 3 This motivated us to construct 
wavelets belo nging to each equivalence class M n After we constructed these wavelets we 
became aware of the article [SW], in which S Schaffer and E Weber constructed wavelets 
for the classes M n ,n > 1 by using the method of “operator interpolation” (see [DL]) 
In fact, they did this for other dilation factors as well In this section we will construct 
wavelets belonging to each of the classes M n , n > 0 by a different method Our approach 
will be simpler than that of [SW] m the sense that for each n > 3, we will construct a 
function ip n such that ip n has the required properties to be in M n -i as characterized by 
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Weber Then we will show that ijj n is a wavelet Further, we will construct a famil y of 
wavelets belonging to the class Mq 

Note that, m view of the characterization of M n , to prove that rp n e M n - 2 , it is 
sufficient to show that ip n satisfies the following two conditions 

(a) ( supp fa) fl ( supp + 2 k qn) = 0, for all q € 2Z + 1 and k = 1, 2, n - 2 

(b) There exist a subset E of supp ip n such that E+2 n ~ 1 qn c supp i> n for some q e 2Z+1 


3.3.1 Construction of wavelets in M n , n > 1 

Let n > 3 and define a n , b n , Cn, d n and e n as in Chapter 2 That is, 

l r\ 2 n 

0*71 2 n — 1^* J — 2 n — 1 ’ 




dn = 


aiii* 
2 n — 1 " » 


e - 

c n — 2 n — 1 ' 

Our starting point is the MSF wavelet rj given via its Fourier transform fj = xw n , where 


W n — ~b n , ~~ o-n. 0 [a n , e n U [c n , d n 


The wavelet 77 belongs to the class of wavelets characterized in Theorem 2 1 of Chapter 2 
(see (B) of §2 6) 

We will construct ip n from this function m the following manner We translate 
(If > if] + 2” _1 7r C [Cn, d n ] to the left by a factor of 2 n ~ 1 n and assign values to ij) n on 
both these sets Then, we translate [a n , e n ] to the right by a factor of 2 n 7r and assign to 
4> n the value ^ on [a n , e n ] and on [a n , e n ] + 2 n n We leave the remaining sets of W n 
as it is More precisely, we have the following function 


1 


^n(£) = < 


1 

1 

“71 


if £ e [-bn, -On] u [Cn,f + 2 n “ 1 7r] U [f + 2 n - 1 7T, d n ] 
if £ € I s ?, f] U [c a n , e„] U [<f + 2 n_1 7r, f + 2 Tl “ 1 7r] 
if £ € [a n + 2 n 7T, e„ + 2 n 7r] 


v 


0 


otherwise 
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Let us make some observations about the set F n = supp 4> n It is easy to verify the 
following facts (stated in a tabular form for easy reference) 


£ € the set 

£ 4- 2/c7r £ P n unless 

2 J £ ^ F n unless 


k = 0 

bbh 


k = Q 

EBBBI 

f 9hl &l1 

L 2 J 2 J 

k = 0, 2 n ~ 2 

J = 0,1 

[ttn? &n\ 

k = 0, 2 n ~ 1 

J = 0,-1 

[Cn, ^ + 2”“%] 

k = 0 

J = 0 

[f + 2 n-1 7r, f + 2 n -V] 

k = 0, — 2 n-2 

J = 0,1 

[f + 2”- 1 7r,d n ] 

/c = 0 

J = 0 

[a n + 2 n 7r, e n + 2 n 7r] 

Jfe = 0, — 2 n-1 

J = 0,-1 


Theorem 3.3. For each n > 3, the function ip n defined above is a wavelet and belongs to 
the equivalence class M . n ~2 

Proof: To prove that ip n is a wavelet, it is sufficient to show that ip satisfies the following 
three properties (see Theorem 1.1) 

(1) HiMIs = i 

(2) ff(() = £ li.(2’?)l 2 = 1 for a e £ e R 

jez 

(3) tq(Z) = £ i’ni 2 , O&( 2, (£ + 2 «* 0 ) =0 for a e £ € R and for all q e 2Z + 1 

J>0 

Proof of (1) 

We have 


lllh.ll? = (6„-o») + (y + 2"-V-o„)+4-(| + 2’- 1 x) 

+ 5{(T-T) + ( ^- aJ+( T-T) + ( ^-^} 
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— b n - a n 4- d n - Cn + e n - o n = 2ix 

Therefore, \\ 1 pnW 2 = 1 
Proof of (2) 

Since H(£f) = , it is enough to show that H(f) — 1 for a e £ E [a, 2a] for some 

a > 0 We will prove that H(£) = 1 for a e £ G [a», 2o n ] = [a n , b n ] = [a n , e„] U [e n , b n ] 
Suppose £ G [o n , e n ] Then 2^£ € F n only when j = 0,-1 So, if(£) = |^(£)| 2 + 

I^©I 2 = (A ) 2 + (A ) 2 = 1 Now - 

f e [e„,6„] <i^2 n_1 £ G [Cn,4] 

= [cn, y + 2 n-1 7r] U [y + 2 n - 1 7T, y + 2""^ 

u[|+2"-Va] 

= J1U/2U/3, say 

If 2 n-1 £ G (/1 U 13), then 2 J (2 n-1 £) € only when j — 0 (from the table) So, H(g) = 
|4(2 n-1 £)| 2 = 1 Also if 2 n_1 £ G I 2 , then 2^(2 n " 1 £) G only when j = 0 or 1 So 
H(0 = |^n(2"- 1 e)| 2 + |^n(2-e)| 2 = (^) 2 + (-^) 2 = 1 Hence, ff(£) = 1 for a e £ > 0 
For £ < 0, it suffices to show that i?(£) = 1 on -a„] In fact, from the table it 
is clear that if £ G [-6 n , then 2 J £ G F n only when j = 0 So, #(£) = |^ n (£)| 2 = 1 
for a e £ G [— b n , — a n ] Thus, we have proved that #(£) = 1 for a e £ G M 
Proof of (3) 

We now have to prove that 

t ? (£) = T/> n (2 J £)i/> n (2J(£ + 2qir)) = 0 for a e £ G R and for all q G 2Z + 1 

j>0 

Since t,(£) = i_ g (£ + 2g7r), it is enough to show that t q = 0 a e , if q is a negative odd 
integer Suppose q —1, and odd We have 2°q ^ 0, ±2 n_1 ± 2 n ~ 2 Therefore, if 2 J £-F n , 
then from the table it is clear that 2 J £ + 2 2 3 qir & F n , which shows that each term of the 
sum t 9 (£) is 0 Hence, t q = 0 a e 

It remains to prove that t_i(£) = 0 for a e. £ € R 

t_i(£) = 2Jjr) 

J>0 
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From the table (see the 8th and 6th row), we observe that both 2 J £ and 2 J £ - 2 2% 
belong to F n only in the following two cases 

(i) 2 J £ G [a„ + 2 n 7r, e n 4- 2 n 7r] and j = n - 1 
(a) 2 J £ G 4- 2 n-1 7r, 4- 2”“%] and j = n — 2 

But both are equivalent to saying that 2 n “ 1 £ e [o n + 2 n ?r, e n 4- 2 n 7r] So we get t_ : (£) = 0, 
if2 n ” 1 £g [o n + 2 n 7r,e n 4-2 n 7r] 

Now if 2 n-1 £ € [a n + 2 n 7r, e n + 2 n 7r], then 

MO = &(2 n ~ a £)$„(2 n " a £ - 2"-%) + ^ n (2 n - 1 £)V-n(2 n - 1 £ - 2”7t) 

= JL J_ 

^ x/2 \ V2j V2 

= 0 

We have proved that t 9 (£) = 0 for a e £ G R and for all q G 2Z 4- 1 Therefore, i/> n is a 
wavelet 

Our claim now is that ip n G M n - 1 For this, we have to show that 

(a) F n fl ( F n 4- 2 k qn) = 0 for all q G 2Z 4- 1 and k = 1,2, n - 2, and 

(b) there exists a subset E of JP n such that E 4- 2 n ~ 1 qir c F n for some q G 2Z 4- 1 

As 2 fc 57 r = 2 2 fc-1 <77r, by referring to the table we observe that if £ G F n , then £4-2 2 k ~ l qir £ 
F n (since q is odd and k — 1 < n - 3) So (a) is proved For (b) consider E = [ 9 f, § ^] and 
q = 1 Then, E 4- 2 n ~ 1 qir = [3f, 4- 2 n_1 7r C F n 

Since n > 3, we have proved that the equivalence classes M n , n > 1 are non-empty 


3.3.2 A Family of Wavelets belonging to the class Mo 

For n > 3, let a n , b n , c n , d n , and e„ be as above In Theorem 2 1 of Chapter 2, we define 
b on [e n , b n ] as follows 


m = 


i 

72 

0 


if £e [e n ,7r] 
if £ G [tt, b n \ 
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Then we extend b to the whole of S n by using (*) - (*«) of Theorem 2 1 Observe that 
[e n , b„\ Pi supp b D (^=r supp b) = [e n ,7r] We define 9 on [e n> ?r] as 

* if £ € [e n ,7r] 

0 otherwise 

Clearly, 6 satisfies the functional equation m (v) of Theorem 2 1 This choice of b and 9 
will give us the wavelet w n , where 

1 if £ G [-7T, -a n ] U [a„, e n ] U [2 n_1 7r, d n ] 

j- 2 if £ € [~d n , -2»- 1 7T] U [-bn, 7r] U [c„, 2"- 1 7T] 

~72 lf ^ 6 M 
0 otherwise 

\ 

Weber [Web, Theorem 4 ] also proved that a wavelet i/i e L\ only if supp i> is not partially 
self similar with respect to any odd multiple of 2w In view of this result, to show that 

A 

w n e Mo = Co \ jCi, we only have to find a set H n m supp ip n such that H n + 2qn is also 
a subset of supp ip n , for some odd integer q The choice H n = [— b n , — tt] and 9=1 will 
do the job, as [- b n , -7r] + 2n = [ e n , 7r] C supp w n 




Chapter 4 


Interval Wavelet Sets of if 2 (R) 

4.1 Introduction 

The classical Hardy space H 2 ( R) is the collection of all functions of L 2 ( R) whose Fourier 
transform is supported in 1+ = (0, oo) 

H 2 ( R) = {/ e T 2 (R) /(£) = 0 for a e £ < 0} 

It is clear that H 2 (R ) is a closed subspace of L 2 (R) As m the case of L 2 (R), we can 
define a wavelet for H 2 (R) A function ip Q. H 2 (R) is sold to be a wavelet of H (R) if the 
system of functions 

{ip h k = 2 j/2 (2 j -fc) 

/onus an orthonormal basis for H 2 { R). Such a function ip will be called an H 2 -wavelet 
An example of an H 2 - wavelet is the function whose Fourier transform is the character- 
istic function of the interval [27T, 4 tt] In fact, for a long time this was the only known 
tf 2 -wavelet P Auscher [Aus] proved that if ip E H 2 { R) is such that \ip\ is continuous 
and | , t/>(£)| = 0((1 + |£| _a ~s) at oo, for some a > 0, then ip cannot be an i7 2 -wavelet In 
particular, there is no band-limited if 2 - wavelet such that \ip\ is continuous 

Similar to the case of L 2 { R), all tf 2 -wavelets are characterized by the following four 


equations 
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^|^( 2 ^)| 2 = Xr + ( 0 for a e £ > 0 (41) 


X ( 2J (f + 2 ?tt)) =0 for a e f e R and for all g e 2Z + 1 (4 2) 

j>0 

Xl^ + 2^1 2 ^ 1 for a e £ € R (4 3) 

A-EZ 


X + 2fc7r )^ ( 2J (£ + 2 M) = 0 for a e £ € R and for all j > 1 (4 4) 

AlEZ 


An 77 2 -wavelet 0 will be called an MSF (minimally supported frequency) wavelet if 
|0| = Xk for some measurable subset K of R+ The associated set K, which has measure 
27 t, will be called an i7 2 -wavelet set It can be proved that (see Theorem 1 2, Chapter 
1 for the L 2 case) a function 0 e H 2 (R) such that |0| = xk for some A c 1+ is an 
JT 2 -wavelet if and only if (4 1) and (4 3) are satisfied We can rewrite the above result in 
terms of the wavelet set 


Theorem 4.1. A set K C1R+ is an H 2 -wavelet set if and only if 


(I) {K + 2&7T k £ Z} is an a e partition of R 
(II) {2 JR j £Z } is an a e partition of R+ 


An H 2 - wavelet set K is said to be an interval H 2 -wavelet set if it is the union of a finite 
number of intervals of R+ In [HKLS] the authors characterized all J7 2 -wavelet sets which 
is either an interval or a union of two intervals In fact, the only if 2 -wavelet set which is 
an interval is [27r, 47 t] The ^-wavelet sets which are union of two disjoint intervals are 
the following sets 

2 r+1 k 2 r+2 (k + 1) 

rrr -nr 

2 r - 1 ’ 2 r+1 - 1 J ’ 


K = 


2(k + 1) 2 k 

2 r + 1 - l 7r ’ 2 r - 1 


-7T 


u 


where r > 0, 0 < k < 2(2 r - 1), r,k£ Z (see Theorem 5 2 of [HKLS]) 

In this chapter we prove a result on the structure of interval wavelet sets of H (R) 
and then characterize all 3-mterval j? 2 - wavelet sets 
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4.2 Interval iJ 2 - wavelet sets 

We have the following theorem regarding interval # 2 -wavelet sets 

Theorem 4.2. Let I s = [o„ b s ], 0 < s < m and K = I 0 U h U • U I m , with 0 < a 0 < 
bo < ax < bi < < a m < b m Then K is an H 2 -wavelet set if and only if 

(1) I 0 U (4(1) - 2km) U (7,(2) - 2/c 2 7t) U U (4( m ) - 2fc m 7r) = [a 0 , o 0 + 2 tt], 

for some permutation a G S m , k x G Z+, 1 < i <m such that the right end point of 
each interval is equal to the left end point of the next interval, and 

(2) I 0 U (2 ri / p( i)) U (2 r2 / p(2) ) U • • U (2 r ”*/p (m) ) = [o 0 , 2a 0 ] 3 

for some permutation p G S m , r, G Z + , 1 < i < m such that the right end point of 
each interval is equal to the left end point of the next interval, 

where S m is the permutation group of { 1,2, • ,m} 

Proof: Suppose K is as m the hypothesis and (1) and (2) hold Then, 

m 

u<* + 2fc7r) = UIK* + 2kir) 

kez fee Zi=o 

m 

= u[{ / o u UM’M- 2 W}+ 2 H 

AigZ *=1 

= | [do; Oo + 2 tt] + 2H | 

k gz k 

= R 

Also, {K+2kn k G Z} is pairwise disjoint Similarly, U je z(2 J iiQ = R+ and {2 °K j G Z} 
is pairwise disjoint Therefore, K is a wavelet set for H 2 (R) by Theorem 4 1 

Conversely, suppose that K = Zo U U I m is an H 2 - wavelet set, where I r = [a r ,b r \, 
0 < r < m So, (I) and (II) of Theorem 4 1 are satisfied Note that for each r, 1 < r < m, 
there is a uni que 4 G Z such that a r + 2k r r G [ao, fl o H - 2 tt] But, a r + 2k r x [floi ^o] 
Because in that case, I 0 n(I r +2km) will contain the interval [a r +2k r ir,b 0 ] if 6 0 < b T +2k r -K, 
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or the interval [a r + 2k r n, b T 4- 2k T n\ if b T + 2k T Tr < bo This will be a contradiction to (I) 
So, a r + 2k r ir G [6 0 , a 0 + 27 t] Note that k r > 0 as a r > b 0 
Claim: I T + 2fc r 7r = [a r , b r ] + 2k r n c [b 0 , a 0 + 27r] 

If not, then bo < a r + 2k T n < ao + 2w < b T + 2A r 7r Therefore, a T + 2 (k r — l) 7 r < 
a,Q < b r + 2 (k r — l)7r But, this implies that Iq fi ( I r + 2{k r — l) 7 r) contains the interval 
[ao, b T + 2 (k r — l)7r] if b r + 2{k r — l)7r < bo, or the interval [ao, bo] if bo < b r + 2(k T — l) 7 r 
In either case, (I) is violated Hence, the claim is proved Now, 

m 

U(jr + 2fc7r) = |j[|J(/'' + 2 fcr )' 

fcez kez r=o 

r ( 4 5 ) 

= U [{ (^° u u u (Im + 2^ m 7r) } + 2k'K 

fee z 

If there is no r, 1 < r < m, such that b Q = a r + 2k r n, then {(Jj. + 2 &i7t) U U (I m -t- 
2k m Tr) + 2kn} will be properly contained in [b 0 , a 0 + 27r] In fact, the set [a 0 , a 0 4- 27 t] \ 
{ (lo U (I i + 2fci7r) U U ( I m + 2k m n)) + 2A:tt} will have positive measure which, in turn, 
will show that (4 5) cannot be an a e. partition of R So, there is an index r, 1 < r < m 
such that bo = a r + 2k r ir Further such an index is unique For, if there exist indices r and 
s with 1 < r, s < m such that b 0 = a r + 2k T n = a s + 2k s n, then ( I r + 2ft r 7r) n (I s + 2k s ir) 
will contain an interval which will again contradict (I) 

Therefore, there is exactly one index i\ such that bo — a n + 2kirc, 1 < ii < m and 
ke TL Now, I 0 U (I tl + 2fci7r) = [a 0 , 6i + 2fci7r] C [a 0 , a 0 + 2n) In a similar manner we can 
show that (the role of bo is now taken by b% + 2kiir) there is exactly one index such that 

6 n + 2fci7r = a l2 + 2&27T, 1 < %2 < m, 12 ^ *1, € Z + 

Similarly, 

b t2 + 2A: 2 7 t = a l3 + 2fc 3 7r, 1 < *3 < m, *3 # *i> *2. h € Z+ 

^ m _i + 2fc m _!7r = a lm + 2k m 7t , 1 < ! m5 £ l i> > ®m- ii k m G Z + 
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Now, b lm + 2/c m 7r has to coincide with a 0 + 2 i r Otherwise [a 0 , a 0 + 27r] \ { (I 0 U (J u + 2fci7r) U 
U (I %m + 2/c m 7r)) 4- 2/c7r} will have positive measure which will contradict (I) So, we 
have proved (1) of the theorem Similarly, by considering dilations by powers of 2 of the 
intervals I r and making use of the partition (II) of M+ , we can prove (2) □ 

4.3 Construction of 3-interval H 2 - wavelet sets 

We want to characterize all wavelet sets K of H 2 (R) consisting of 3 intervals In view of 
Theorem 4 2, we have the following result for such a K 

Corollary 4 . 1 . Let K = [ax, &i] U [02, 62] U [03, 63] such that 0 < ai < 61 < a 2 < b 2 < a$ < 
63 Then, K is an H 2 -wavelet set if and only if for some non-negative integers r, s, k, l 

( 1 ) either [Tl] [oi, &i] U ([02, 62] - 2 far) U ([a 3 , 63] - 2 hr) = [01, 01 + 27 r] 

with bi = a2 — 2&7T, b 2 — 2 kir = 03 — 21k, 63 — 2 lir = 01 + 27 r 
or [T2] [oi, 61] U ([a 3 , 63] - 27 d) U ([o 2 , b 2 ] - 27rk) = [a x , 01 + 27 r] 

with bi = a$ — 2/7T, 63 — 2 lr = a 2 — 2 kir, b 2 — 2 kir = Oi 4 - 2 ir 
and ( 2 ) either [Dl] [ai, 61] U 2" r [a2, b 2 ] U 2 _i [a3, 63] = [oj, 2 ai] 

with bi — 2 ~ r a 2 , 2 ~ r b 2 = 2 _s a 3 , 2 _s 6 3 = 2ai 
or [D 2 ] K, U 2- s [a 3 , 63 ] U 2- r [o 2 , b 2 ] = [oi, 2aJ 
with b\ = 2 ~ a as, 2 ~ 3 h — 2 ~ r a 2 , 2 ~ r b 2 = 2ai 

Thus, in order to characterize all H 2 -w avelet sets consisting of three intervals, we have 
to consider each of the four cases (T t ,D,), t,J = 1,2 and determine the values of the 
non-negative integers r,s,k,l such that the corresponding relations (T t , Df) hold 

THE CASE (T1,D1) 

We have K — [ai, 61 ] U [a 2 , b 2 ] U [ 03 , 63 ] 


b\ — U 2 2krr , 

bi = 2~ r a 2 , 


b 2 — 2kir = 03 — 2I7T, 
2 r b 2 — 2 ®3) 


63 — 21tt = Oi + 27T 
2 ~ s b z = 2ai 


(4 6) 
(4 7) 
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Since bi < as, we have ft > 1 Now, b 2 ~ 2kn = a 3 - 21 tt =>a 3 -b 2 = 2(1 - k)ir Since 

a 3 -b 2 > 0 , we have l > k Thus, 1 > k > 1 Similarly, s > r > 1 Solving (4 6) and (4 7) 
for a u b u we get 


“1 = 

h = jSlTT, 

02 = 


“3 = 

= 


We have to ensure, 0 < a x < b x < a 2 < b 2 < a 3 < b 3 Clearly, a x > 0, b x < a 2 and b 2 < 
The conditions a x < b x , a 2 < b 2 , and a 3 < 63 will imply 




( 1 ) (2 r — 1)1 < (2 S+1 — l)ft — (2 r - 1), 


(II) ( 2 r - 1 )/ > ( 2 s - 1 ) k , 

(III) ( 2 r+1 - 1 )/ < ( 2 3+l - l)k + 2 ( 2 s - 2 r ) 

Also, we can show that (111) =>(1) So, we have to consider only (11) and (111) Eliminating 
l from (11) and (111), we get 1 < k < 2(2 r - 1) 

Thus, to get all if 2 -wavelet sets in this case, we proceed as follows 
Fix r > 1, then we have to consider only those k such that 1 < k < 2(2 r — 1) Consider 
any such k By fixing an s > r + 1, we determine all integer l satisfying (11) and (111) 
Then, any such combination of r, ft, s, l will give rise to a wavelet set 

For example, let r = 1, then 1 < k < 2 So, k = 1 If s = 2, then (11) and (111) give 
us l > 3 and 31 < 11 This implies l > 4 and l < 3 So, there is no l satisfying (11) and 
(111) If s = 3, then (11) and (111) imply l > 7 and 31 < 27 =*► l > 8, 1 < 8 =r> l = 8 So, 
r = 1, ft = 1, s = 3, 1 = 8 give rise to a wavelet set. The corresponding wavelet set is 


[7 r + |7r, 27r] U [47 t, 47 r + |7r] U [187T + |7r, 197 t + |7r] 


If 5 = 4, then we get, 1 > 15 and 31 < 59 =>• 1 > 16, 1 < 19 =>• 1 — 16, 17, 18, 19 So if 
we take r = 1 , ft = 1 , s = 4, we get wavelet sets for each of the 1’s, 1 = 16, 17, 18, 19 The 
wavelet set corresponding to r = 1, ft = 1, s = 4, 1 = 16 is 


[tt + frr, 2tt] U [4tt, 4t r + frr] U [34tt + fir, 35tt + £tt] ^ 

WTTcfbT 5pT^| 
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A short table of (r, k,l, s) is given in Table 1 at the end of this chapter 

THE CASE (T2,D2) 

We have K = [a x , 6^ U [a 2 , 6 2 ] U [a 3) b 3 ] and 

bi=a 3 -2l7r, b 3 - 2lir = o 2 — 2kw, b 2 ~ 2kir = a x + 2tt, (4 8) 

b x = 2~ s a 3 , 2~ s b 3 = 2~ r a 2 , 2~ r b 2 = 2 a a (4 9) 


Examining the relations among o t ’s and 6 t ’s as m the case (Tl.Dl), we have l > k > 0 
and s > r > 0 Solving the equations (4 8) and (4 9) for a,’s and 6,’s , we get 


a i — 2 r+i"^i 7r i h = 
a 2 = 

a 3 = m IT, 63 = 

Again, we have to ensure that 0 < a x < b x < a 2 < b 2 < a 3 < b 3 Clearly, 0 < a x As 
r + 1 < s, ai < b x =$■ b 2 < a 3 Also, a 3 < b 3 =$■ b x < a 2 So, we have to consider the 
inequalities a x < b x ,a 2 < b 2 ,a 3 < b 3 These conditions are equivalent to the following 
inequalities 


( I ) (2 r+1 -l)t>(2 s -l)(k + l) 

( II ) ( 2 r+1 - 1)1 < ( 2* +1 - 1 )k + 2(2 S - 2 r ) 

( III ) (2 r - 1)1 > (2 s - 1)A: 

If k = 0, then ( 111 ) is trivially satisfied So, we have to consider only ( 1 ) and ( 11 ) And, 
if k > 0, then ( 111 )=* ( 1 ) So, one has to consider ( 11 ) and ( 111 ) Also, (u) and ( 111 ) imply 
k < 2(2 r - 1) If r = 0, then k < 0, which is not possible. Thus, to get all wavelet sets in 
this case, we proceed as follows 

Fix r > 1 and consider all k such that 0 < k < 2(2 r — 1) Take s > r If k = 0, 
determine all l satisfying ( 1 ) and ( 11 ), and if k > 0, then determine all l which satisfy ( 11 ) 
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and (in) For example, r — 1, /c — 0, s ~ 2 does not give any wavelet set But if we take 
r = l,k — 0,s = 3, then we get l = 3 The corresponding wavelet set is 

[| 7r > H U [2 tt, 2tt + |tt] U [6tt + |tt, 8tt] 

Observe that when k > 0, the inequalities to be considered are same as in the case 
(T1,D1) So, the table for (T1,D1) also works for (T2,D2) though we will get different 
wavelet sets A short table for the case fc = 0 is given in Table 2 

THE CASE (T2,D1) 

In this case, we have, s > r > 0, l > k > 0. Solving the equations [T2] and [Dl] of 
Corollary 1 for a t ’s and 6,’s, we get 

Oi = 2 ^t[( 2 3 — ~ (^ r — + 2 s ]7T , 

b x = 2^t[(2 6+1 — l)k — (2 r+1 — 1)1 + 2 s+1 ]7r, 

a 2 = 2[(2 a+1 — 1)A; — (2 r+1 — 1)1 + 2 a+1 ]7r, 

b 2 = 2^r[(2 s+1 -l)fe-( 2 r -l)M- 2 s+ 1 ] 7 r, 
a 3 = ^[( 2 s+ 1 -l)fc-( 2 r -l)I + 2 s+ >, 

6 3 = 4[(2 s -l)fc-(2 r -l)/ + 2 3 ]7T 

Clearly, b x < a 2 and b 2 < 03 The conditions 0 < a x ,a x < b x , a 2 < b 2 and 03 < 63 will 
imply 

(I) (2 4 - 1 )k + 2 S > (2 r - 1)1 

(II) [2 s (2 i+1 - 1) - 2 r (2 s - 1 )]k + (2 23+1 - 2 r+a ) > [2®(2 r+1 - 1) - 2 r (2* - 1)]/ 

(III) (2 s - 1)(2 S+1 - l)fc + (2 2s+1 - 2 s - 1 -. 1 ) < [2 s (2 r+1 - 1) - (2 r - 1)]/ 

(iv) [2 r+1 (2 a - 1) - (2 S+1 -l)]k + 2* +1 (2 r - 1) > (2 r - l)(2 r+1 - 1)1 

One can show (11) =* (iv) => (l) So, we have to consider only (11) and (111) Eliminating l 
from (111) and (iv), we get k < 2(2 r - 1) To get all wavelet sets in this case, we apply the 
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similar procedure adopted in the case (T1,D1) See Table 3 for some acceptable values of 
(r, k, s , 0 

THE CASE (T1,D2) 

Here we have, s > r > 0, l > k "> 1 Solving [Tl] and [D2] for a^s and 6 t ’s, we get 

a i = [(^ r — ' — (2 s — 1)& + 2 r ]7T, 

h = —=x [(2 r+1 — 1)1 — (2 S+1 — l)k + 2 r+1 ]7r, 

°2 = ^[(2 r+1 -l)Z-(2*-l)fe + 2 r+1 ]7r, 

b 2 = 4[(2 r -l)l-(2 s -l)k + 2 r ]w, 
o 3 = 2[(2 r+1 — 1)1 — (2 i+1 — l)k + 2 r+1 ] 7r, 
h = ^[(2 r+1 -l)l-(2 3 -l)k + 2 r+1 ]n 

The conditions 0 < Ox, ox < &x, a 2 < b 2 and o 3 < & 3 are equivalent to the following 
inequalities 

(l) (2 r - 1)1 + 2 T > (2 s - l)k 

(n) [2 s (2 r - 1) - 2 r (2 r+1 - 1 )]l + ( 2 r+s - 2 2r+1 < [2 s (2 s - 1) - 2 r (2 s+1 - 1 )]k 

(m) [2 s+1 (2 r - 1) - (2 r+1 - l)]l + 2 r+1 (2 s - 1) > ( 2 s - 1)(2 S+1 - l)fc 

(iv) (2 r - l)(2 r+1 - 1)1 + 2 r+1 (2 r - 1) < [2 r (2 4+1 - 1) - (2 a - 1)] k 

The coefficient of l in (u) is negative if and only if s = r + 1, in which case (u) is trivially 
satisfied If this coefficient is non-negative, then it can be shown that (iv)=» (u) Also, it 
can be shown that (m)=» (i) So, we need only (m) and (iv) As (m) and (iv) imply that 

k < 2(2 r — 1), the case r = 0 is ruled out because k > 1 Table 4 gives few acceptable 

values of (r, k, s, l ) 
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Table 1 

Table 2 

Table 3 

Table 4 

(T1.D1) 

(T2,D2), k = 0 

(T2,D1) 

(T1,D2; 


r 

k 

s 

l 


—1 





6 


5 3,4 


3 8 

4 16-19 

5 32-40 

6 64-83 
128-168 
256-339 
512-680 

1024-1363 


3 

6,7 

11-16 

22-35 

43-71 


5 

11,12 

21-25 

43-53 
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11 3 

12 5-7 

20 1026-2047 

25 32801-65567 


33-63 


1 2 
3 


3 3 

4 6-9 

5 11-19 

6 22-41 

7 43-83 


2 3 
4 



5 6 
7 


7 3 

8 5-7 

9 9-15 

10 17-31 











































Chapter 5 


Multiwavelet packets and frame 
packets of L 2 (R d ) 


5.1 Introduction 

Consider an orthonormal wavelet of L 2 (R) At the j-th. resolution level, the orthonormal 
basis generated by the wavelet has a frequency localization proportional to 2 J For exam- 
ple, if ^ is a band-limited wavelet, then the measure of the support of ('ip J k) A is 2 J times 

A 

the measure of the support of since 

(M A (0 = 2-*$(r'0e- 2_, «, 3 , k e z 

So when j is large, the wavelet bases have poor frequency localization Better frequency 
localization can be achieved by a suitable construction starting from an MRA wavelet 
basis 

Let {V 3 j € Z} be an MRA of L 2 ( R) with corresponding scaling function <p and 
wavelet ijj Let W } be the corresponding wavelet subspaces W 0 = sp{ipjk k G Z} In 
the construction of a wavelet from an MRA, essentially the space Vi was split into two 
orthogonal components Vo and W 0 Note that Vi is the closure of the linear span of the 
functions {2 s</j(2 - k ) k G Z}, whereas V 0 and W 0 are respectively the closure of the 
span of -{X - k) k} and {V>( - k ) k} Since <p( 2 - k ) = ip( 2( - |)), we see that the 
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above procedure splits the half integer translates of a function into integer translates of 
two functions 

In fact, the splitting is not confined to V[ alone we can choose to split Wj, which is 
the span of (ip( 2 3 — k ) k} — {ip (%>( — &}, to get two functions whose 2 

translates will span the same space Wj Repeating the splitting procedure j tim es, we 
get 2 ] functions whose integer translates alone span the space W 3 If we apply this to 
each Wj, then the resulting basis of L 2 (R), which will consist of integer translates of a 
countable number of functions (instead of all dilations and translations of the wavelet ip), 
will give us a better frequency localization This basis is called “wavelet packet basis” 
The concept of wavelet packets was introduced by Coifman, Meyer and Wickerhauser 
[CMW1, CMW2] For a nice exposition of wavelet packets of L 2 ( R) with dilation 2, see 
[HW] 

The concept of wavelet packet was subsequently generalized to R d by taking tensor 
products [CM] The non-tensor product version is due to Shen [She] Other notable gen- 
eralizations are the biorthogonal wavelet packets [CD], non-orthogonal version of wavelet 
packets [CL], the wavelet frame packets [Che] on R for dilation 2, and the orthogonal, 
biorthogonal and frame packets on R d by Long and Chen [LC] for the dyadic dilation 

In this chapter we generalize these concepts to R d for arbitrary dilation matrices and 
we will not restrict ourselves to one scaling function We will consider the case of those 
MRAs for which the central space is generated by several scaling functions 

Definition 5.1. A dx d matrix A is said to be a dilation matrix for R d if 

(i) A(Z d ) C Z d and 

(it) all eigenvalues X of A satisfy |A| > 1 

Property (i) implies that A has integer entries and hence | det A\ is an integer, and (m) 
says that the integer is greater than 1 Let B = A t , the transpose of A and a = | det A\ = 
| det JB 1 
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Considering Z d as an additive group, we see that AZ d is a normal subgroup of Z d So 
we can form the cosets of AZ d in Z d It is a well known fact that the number of distinct 
cosets of AZ d in Z d is equal to a = | detA| ([GM], [Woj]) A subset of Z d which consists of 
exactly one element from each of the a cosets of AZ d in Z d will be called a set of digits 
for the dilation matrix A Therefore, if Ka is a set of digits for A, then we can write 

Z d = |J (AZ d + p), 

where {AZ d + p p £ Ka} are pairwise disjoint A set of digits for A need not be a set 
of digits for its transpose For example, for the dilation matrix M = (Jo) of R 2 , the set 
{ (o) ’ (o) } 1S a set °f digits for M but not for M* It is easy to see that if A is a set of 
digits for A , then so is K - p, where p £ K Therefore, we can assume, without loss of 
generality, that 0 £ K 

The notion of a multiresolution analysis can be extended to L 2 (R d ) by replacing the 
dyadic dilation by a dilation matrix and allowing the resolution spaces to be spanned by 
more than one scaling function 

Definition 5.2. A sequence { V 3 . j £ Z} of closed subspaces of L 2 (R d ) will be called a 
multiresolution analysis (MRA) of L 2 ( R d ) of multiplicity L associated with the dilation 
matrix A if the following conditions are satisfied 

(Ml) V, C V 3+1 for all j £ Z 

(M2) Ujg zVj is dense m L 2 (R d ) and C\ j€ zVj = {0} 

(MS) f £ V 3 if and only if f (A) £ V J+ i 

(M4) there exist L functions {^i, </? 2 , . <Pl) m v o, called the waling functions, such that 

— k) 1 < l < L, k £ Z d } is an orthonormal basis of Vo 

The concept of multiplicity was introduced by Herve [Her] m his Ph D thesis Since 
1 < l < L, k £ Z d } is an orthonormal basis of Vo, it follows from property 
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(M3) that {a 3 / k) , 1 < / < L, k G Z d } is an orthonormal basis of V, Observe 

that if / € L 2 (R d ), then 

(a j/2 f(A -A;)) a (0 = o -J / 2 e - *( s "^ ,fe )/(5 - i^), ( E R d , k E Z d 

First of all we will prove a lemma, the splitting lemma (see [Dau]), which is essential 
for the construction of wavelet packets We need the following facts for the proof of the 
splitting lemma 

(1) Let T* = [-7r,7r] d and / € L 1 ^) Since M d = U JfeeZ < i (T d + 2 tor), we can write 

f f(x)dx = [ { X) f( x + 2kn)\dx (51) 

J J T d '■fcg Id i 

(2) Let {sk k E Z d } E l x ( Z d ) and Kb be a set of digits for the dilation matrix B As 
Z d can be decomposed as Z d = U ^k b (BZ d + p ) , we can write 

E s k = E E s H+Bk (5 2) 

fc€ z d I*€Kb kez d 

(3) Let Kb be a set of digits for B Define 

Qo= |J B~ 1 (T d + 2p%) 

)i€Kb 

Since Kb is a set of digits for B, the set Qq satisfies ^Jkez d {Qo + 2kn) = M d This 
fact, together with |Qo| = 1> implies that {Qo + 2&7 r k E Z d } is a pairwise disjoint 
collection (see Lemma 1 of [GM]) Therefore, 

[ f(x)dx = [ \j:f{x + 2k7r)}dx, for / G L 1 (R d ) (5 3) 

J R d JQo '•/fee z d J 

A function / is said to be 2?rZ d - periodic if f(x + 2kir) = f(x) for all k E Z d and for 
ae £ 

5.2 The splitting lemma 

Let {ipi 1 < l < L} be functions in L 2 (R d ) such that {</?j( — k) 1 < l < L, k E Z d } 
is an orthonormal system. Let V = sp{a^ 2 <pi(A —k) l,k} For 1 < 1,3 < L and 
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0 < r < o — 1, suppose that there exist sequences {h\ Jk k € Z d } e l 2 (Z d ) Define 

ff ( x ) = E E h r ljk a 1/2 <Pj(Ax - k) (5 4) 

3 = i fee z d 

Taking Fourier transform of both sides 

fm = E E 

j=ifcez d 

= E (5 5) 

J=1 

where 

AL({) = E 1 < 1, J < i, 0 < r < a - 1, (5 6) 

fee z d 

and ^ is 27rZ d -penodic and is in L 2 ^) Now, for 0 < r < a - 1, define the L x L 


matrices 



(5 7) 

By denoting 

$(a) = »¥>£(*))* 

(5 8) 

$(£) = (&(£)> 

(5 9) 

we can write (5 5) as 

^r(e) = -ffr(B- 1 0^ _1 ^), 0 < r < a — 1, 

(510) 


where F r (i) = (JT(*), >?(*)• > «(*))' The foUowmg lemme charactemes the orthonor- 

mality of {v’lt* — A) ' 1 < i < A € % A ) 

Lemma 5.1. The system {*>,(- A) 1 < l < L, k € Z"} « orftonomo! */ only if 

+ 2fc?r ) = ^ 

fcez d 

Proof: The proof is similar to the proof of Proposition 1 1 of Chapter 1 
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Lemma 5.2. (The splitting 

Lei {# 1 — ^ ^ *•} 1,1 J uhnil U " * >l L 2 {R d ) such that the system {a^tp^A - k ) 

1 - J - L ' * 6 ^ 1 let V be its closed linear span Let K be a set of 
digits for B Also let ff, H r h< l{ ' Then 

{fi(~ k) 0 ^ >' < o-l, 1 <!<£,, keZ d } 
is an 07 thononnal system if »»«»/ nnhj 

E MM f “ >/# % Mr) + =8 rs Ii (5 11) 

Moreover, {ff('-k) ij,k} » (lH 0r thonormal basis of V whenever it is orthonormal 
Proof: Foz 1 < l,j < L, 0 i « * fa a ■— l and p € Z d , we have 

</;.//( -p)> 

= (2^<(/;) a .(/,’( -?»'') = (2^3 JjmoWW^is 

= ^(B^MB' 1 ?))^ 0 ^ (by(5 5)) 

= jOT / £ EE {AJ»(S->( f + 2 fc,)®B-‘K + 2fcr)) 

^7r; A.€Z rf w*=l n^l *• 

+ 2t I r))SC g::T (f : i r 2fe))} e'W +21 ' ) ^ (by (5 D) 

= tAj/ EEEE AJmO-'tf + 2(BA + + 2 (Bk + ,»)«•)) 

Jjd^ K AeZ d r „-i n ssi 

+ 2(Bk + /^) 7r ))^n(B~ 1 (^ + 2(PF : r /ijTrjye^df (by (5 2)) 

= p/eee + 2B-»®B =i { + 2B_ >) 

{ E M5 _1 (£ + 2^*0 + 2 ^)^(£-‘ I (iF^ r ) + 2M je l( ^>d£ 

= ?2^W f EEE ^jrn(^ l C+2B-V^)^( jB " 1 ^ +2B ~ lM7r ^ rnnel< ^ > ^ 

(by Lemma 5 1) 

= p/jEE 
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Therefore, 


</;>/*( -P)> = 6 ra 5 3l 6<> p 

^ E E ^jm(^ *£ + 2S 1 /x7r)/if m (.B~ 1 £ + 2J5~ 1 /i7r) = forae £ 6 R d 

m€A m=l 

^ E E ^ r m (£ + 2B- : V)h? m (£ + 2 B- Vtt) = <5 rs fy forae £ € R d 

/iGA' m=l 

^ E ■ffr(£ + 2S~V7i')^(^ + 2 ^~V7r) = 5 rs /£, forae £ € R d 

lx€K 


We have proved the first part of the lemma 

Now assume that {/[( - k) 0 < r < a - 1, 1 < l < L, k e Z d } is an orthonormal 
system We want to show that this is an orthonormal basis of V Let / e V So there 
exists {c JP p € Z d } € l 2 (Z d ), 1 < j < L such that 


/(*) = E E c 3P a 1/2 <p 3 (Ax - p) 

j=lpgZ d 

Assume that f -L f[{ — k) for all r, l, k 
Claim: / = 0 

For all r, l, k such that 0 < r < a — 1, 1 < l < L, k G Z d , we have 

o = </f( - *), /) = (/,'( - *), E E CnfP'vM -p)) 

' J=lpez d 

= 75^3 [ E E 

(27r) 7 R d ;=ipez d 

= ££ f £ WS.-ME E (by (5 5)) 

(2tt) m=l J— lp€% d 

= 7^/ E *S»(f)(MOE E 

(27rj jRd m= x .?=ipez d 

= 7^3/ ee«+*8*i 

(27T) JQ 0 qQZim-1 

E E C^^(£ + 2<?tt) e- l( *’ B(i+297r)) e I < Pl?+2 ^ > d£ (by (5 3)) 

J=lp6Z d 

= 7^3 [ E E E E + 2«7r)}e-<*' B «e > <M> I i5 

(2n) d jQ 0 r£ij=i P ez* Uez d > 
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r ^ 

= FhrY /,, ^ ^m(O c mpe *^> s f) e % Md£ (by Lemma 5 1) 

\**'W J Qo m=l pgjd ' 

o}' 2 r l 

~ ToZy ^ / SI] h\m (0 e l(PiC > 

V"'*/ J fl“ 1 ('T <i +2/iir) m=l pggd 

a l/2 /■ ^ 

= (5?j& Lr< + ^ V^e^-V, 

a 1 / 2 r ( L 

= E^y / , „1 £ £ E ^ m (f+2B- 1 H^e‘(> i 'S+ 2S -v)U-(*.Bs ( if 

J B-'T d KneK m=lpgzd J 


Smce |(|^re 1 ^> B > A: € Z d j is an orthonormal basis for L 2 (J3 1 T d ), the above equations 
give 

EEE ^e , ( f+2s "'’’>’>A[ m K + 2B-» = 0 a e for all r, I 

ms= I pg%d 

Form = 1,2, L, define 

C m ({) = E (5 12) 

pez d 

So we have 

E E C m (£ + 2B- V7r)h[ m (^ + 2 £“ V) = 0, 0<r<a-l, 1<Z<L (5 13) 

txZK m= 1 

Equations (5 11) are equivalent to saying that for 0 < r < a - 1, 1 < l < L and for a e 
(el* 1 , the vectors 

+ 2B~ V) 1 < m < L, p € if) 

are mutually orthogonal and each has norm 1, considered as a vector in the aL-dimensional 
vector space G aL , so that they form an orthonormal basis for C aL Equation (5 13) says 
that the vector 

(c m {£ + 2B~V) 1 < m < L, !X € K) (5 14) 

is orthogonal to each member of the above orthonormal basis of C ai Hence, the vector 
m the expression (5 14) is zero In particular, C m (£) = 0, for all m, 1 <m< L That is, 
Cpjp = 0, 1 < m < L, p e Z d Therefore, / = 0 This ends the proof □ 

The splitting lemma can be used to decompose an arbitrary Hilbert space into mutually 
orthogonal subspaces, as m [CMW2] We will use the following corollary later 
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Corollary 5.1, Let { E A 1 < l < L, k 6 Z d } 6e an orthonormal basis of a Hilbert space 
H Let H r , 0 < t < a — 1 £>e as above and satisfy (5 11) Define 

^ = E EWA' 0<r<a-l, 1<1<L, A; G Z d 

T7i=l pgZ^ 

Then { 1 < / < L, A; € Z d } is an orthonormal basis for its closed linear span TL and 

n = e“io^ r 

Proof: Let <^i , </? 2 , , <Pl be functions m L 2 (R d ) such that {</?*( - k ) 1 < l < L, k e Z d } 
is an orthonormal system Let V = sp{a 1//2 </?i(yl —A) l, A;} Define a linear operator T 
from the Hilbert space V to % by T(a 1 ^ 2 ipi(A —k)) = Ei )k Let f[ are as in (5 4) Then, 
T(f[( — k)) — Ff k Now the corollary follows from the splitting lemma □ 

5.3 Construction of multiwavelet packets 

Let {Vj j 6 Z} be an MRA of L 2 (R d ) of multiplicity L associated with the dilation 
matrix A Let {ipi . 1 < l < L} be the scaling functions Since <pi, 1 < l < L are in 
Vo C Vi and {a l/2 (p 3 (A --k), l<J<L,keZ d } forms an orthonormal basis of Vi, there 
exist {h hk ke Z d } € l 2 (Z d ), for 1 < l ,J < L such that 

<pi(x) = E E hi ]k a 1/2 ipj{Ax - k) 

3= ikez d 

Taking Fourier transform, we get 

wK) = EE 

3= 1 k€ Z d 

= ZhuiB-'OhiB-'t), (5 15) 

1 

where h h (£) = £ a~ l l 2 h ljk e-^' k) , and h 3 is 27rZ d -penodic and is m £ 2 (T*) Let J3o(0 

kez* 

be the L x L matrix defined by 


#o(0 - ((M£)) 1$Ij5;l 
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We will call Hq the low pass filter matrix Rewntting (5 15) in the vector notation 
(5 8) and (5 9), we have 

<£(0 = (5 16) 

Let Wj be the wavelet subspaces, the orthogonal complement of V 3 in V 3+ i 

W } = V 3+l QV 3 

Properties (Ml) and (M3) of Definition 5.2 now imply that 

w,±w f , 

and 

feW 3 &f{A->)eWo (517) 

Moreover, by (M2), L 2 (R d ) can be decomposed into orthogonal direct sums as 

£ 2 (K') = (518) 

jez 

= V 0 ®(©W,) (519) 

J>0 

By Lemma 5 1 and equation (5 15), we have (for 1 < 1,3 < L) 

6,t = E ft(f + 2 

k<Z Z d 

= E {EM jB " I K + 2fer))^(^ 1 (f + 2 * tt))} 

ke.z d l, »n=i 

{E MFHe + 2 kK))ipn{B-\t + 2forj)} 

Now, using (5 2), we have 

S jt = E £ £ Vt 6 "^ + 25-V7r)MF T ? + 2FV) 

l^eKs m=l n=l 

E l^ m (5 -1 (^+ 2/«r) + 2for)f£„(.B~ 1 (£ + 2 /itt) + 2A:tt) j- , 
fcez d *• 

where RTg is a set of digits for B Usmg Lemma 5.1 again, we get 

S 3 i= E £ + 2B- l pit)hi m {B~H + 2 5- V) ( 5 20 ) 

/*€!?£ ™=1 
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This is equivalent to saying that 

£ #o(£ + = I L foraef 

Equation (5 20) is also equivalent to the orthonormality of the vectors 

(/r 0 (£ + 2B" 1 /L«r) 1 < j < L, p G K^j, 1 < l < L, £ e T* 

These L orthonormal vectors m the aL-dimensional space <Cf L can be completed, by 
Gram-Schrmdt orthonormalization process, to produce an orthonormal basis for C aL Let 
us denote the new vectors by 

"b 25 * /itt) 1 5* J ^ L, // CE Kb^ i 1 ^ 5: L, — 1, £ £ T^, 

and extend the functions h\ 3 , 1 < r < a - 1, 1 <l,J <L 27rZ d -periodically (see [GLT]) 
Denoting by # r (f), 1 < r < a - 1 the X x L matnx 

we have 

£ iX r (^ + 2B- 1 H^;(^ + 25-V7r) = 5 rs Ii < for a e £ 

USKb 

Now, for 1 < r < o — 1, 1 < l < L, define 

#«) = < 521 > 

J=1 

Since h r l;) are 27rZ d -periodic, there exist A: £ Z d } £ ^ 2 (Z d ) such that 

kez d 

Now, applying the splitting lemma to VI, we see that {f[ ( — k) 0 < r < a — 1,1 < l < L, 
k € Z d } is an orthonormal basis for V x We use the convention <pi = f°, 1<1<L with 
hij = ftfj and hi ok = h^ k The decomposition V\ = Vo © Wo, and the fact that {/“(• - k) 

1 < i l < L, fee Z d } is an orthonormal basis of Vo, imply that 

{//•(. - k) l<r<a-l,l<l<L,keZ d } 
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is an orthonormal basis for W 0 By (5 17) and (5 18), we see that 

{a> /2 f[(A> - k ) 1 < r < a - 1, 1 < l < L, j e Z, k € Z d ) 

is an orthonoimal basis for L 2 ( R d ) This basis is called the multiwavelet basis and the 
functions {/) 1 — r — Q. ~ 1, 1 < l < L} are the multiwavelets associated with the MRA 

{Vj j € 1} For 0 < r < o- 1, by denoting F r (x) = (/[(a;),/ 2 r (a;), ,fl{x))* and 

= (/i ( 0 . f‘i (0» . we can write (5 16) and (5 21) as 

F r (0 = 0 < r < a - 1 (5 22) 

This equation is known as the scaling relation satisfied by the scaling functions (r = 0) 
and the multiwavelets (l<r<a-l) 

As we observed, by applying the splitting lemma to the space Vi = sp{a 1/2 cpi(A - k ) 

1 < l < L, k € Z d j, we get the functions ff, 0 < r < a — 1, 1 < / < L Now, for any 
n e N 0 = Nu{0}, we define /”, 1 < l < L recursively as follows Suppose that f{, 
r € No, 1 < l < L are defined already Then define 

M = EE AJ*o I/a JT(A I - *), 0 < < < a - 1, 1 < ( < L (5 23) 

J=1*6Z<' 

Taking Fourier transform 

uy+TK) = e (5 24) 

j=i 

In vector notation, (5 24) can be written as 

(f’-H-.r)''© = lUB-'OMB-'O (5 25) 

Note that (5 23) defines /” for every non-negative integer n and every l such that 1 < Z < L 
Observe that /,° = ipi, 1 < l < L are the scaling functions and ff ,1 <r < a-1,1 < l < L 
are the multiwavelets So this definition is consistent with the scaling relation (5 22) 
satisfied by the scaling functions and the multiwavelets 

Definition 5.3. The functions {ff n > 0, 1 < l < L} as defined above will be called the 
basic multiwavelet packets corresponding to the MRA {Vj j € Z} of L (R ) of 
multiplicity L associated with the dilation A 
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The Fourier transforms of the multiwavelet packets 

Our aim is to find an expression for the Fourier transform of the basic multiwavelet 
packets in teims of the Fourier transform of the scaling functions For an integer n > 1, 
we consider the unique “a-adic expansion” (1 e , expansion m the base a) 

7T, = juj H" piU + Q? d* ' H" PjOp j (5 26) 

where 0 < /z, < a - 1 for all i = 1, 2, . J and # 0 

If n can be expressed as in (5 26) then we will say n has a-adic length j We claim 

that if n has length j and has expansion (5 26), then 

F„(0 = H h ( 5 2? ) 

so that is the l - th component of the column vector in the right hand side of 

(5 27) We will prove the claim by induction 

From (5 22) we see that the claim is true for all n of length 1 Assume it for length j 

Then an integer m of a-adic length j + 1 is of the form m — p + an, where 0 < p<a 
and n has length j Suppose n has the expansion (5 26) Then from (5 25) and (5 27), 

we have 


(F m nc = ( f „+„) a k ) 

= h„(b-w„(*- 2 o ^( B ( - )+1) f H (^ f ) 

Smce m = „ + m m + • + ft* *»<« has the de81red fonn * 

induction is complete 

The first theorem regarding the multiwavelet paekets is the Mowing 

Theorem 5.1. Let {/," . n > 0, 1 < l < Q the 
a&ove Then 

« ^ r, k € is an orthonormal basts of 
a? < n < a J+1 — 1? 1 * ^ L ’ * € 1 

Wj» J > 0 
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(n) {fi ( 0 — 71 ^ a 1) 1 < l < L, k e Z d } is an orthonormal basts of 

V,* ) >° 

(m) {/"( “ M ^ °> 1 ^ 1 ^ L > fc ^ Z d } is an orthonormal basis of L 2 [R d ) 

Proof* Since {/" 1 < n < a - 1, 1 < i < L} are the multiwavelets, their Z d -translates 
form an orthonormal basis for W Q So (») is verified for j = 0 Assume for j We will prove 
for j+l By assumption, the functions {/, n ( -fc) a? <n< a J+1 -l, 1 <1<L, k eZ d ) 
is an oithonoi mal basis of W 3 Since / € W 3 f (A-) g the system of functions 

{a 1/2 /”(i4 -&) a J < n < a 3+l - 1, 1 <1<L, ke Z d } 

is an orthonormal basis of Wj+i Let 

E n = spW /2 ff(A--k) 1<1 < L, keZ d } 


Hence, 


W J+ 1 = © E, 

71= OJ 


(5 28) 


Applying the splitting lemma to £„, we get the functions 

9 ,"’ r (x) = E E f^ mt a^C(Ax - k ) ( 0 < r < a - 1, 1 < I < I) (5 29) 

m=l ke z d 

so that {^ n,r ( — k) 0 < r < a - 1, 1 <1 <L, k e Z d } is an orthonormal basis of i? n 
But by (5 23) 

„n,r __ rr+an 

9l ~ Jl 

This fact, together with (5 28), shows that 


{/ H-an ( . _ o < r < o - 1, 1 < l < L, k eZ d , a 3 <n< a 3+1 - 1} 
= {/"( - *) flJ+1 < n < a 3+2 - 1, 1<1<L, ke Z d } 


is an orthonormal basis for W J+ i So (t) is proved Item ( n ) follows from the observation 
that Vj = V 0 © W 0 © • ® Wj-i and (m) follows from the fact that liV 3 = L 2 (R d ) □ 
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5.4 Construction of orthonormal bases from the mul- 
tiwavelet packets 

We now take all dilations by the matrix A and all Z d -translations of the basic multiwavelet 
packet functions 

Definition 5.4. Let {/" n > 0, 1 < l < L} be the baste multiwavelet packets The 
collection of functions 

V = {a j/2 f?{A J - k ) n > 0, 1 < l < L, j E Z, keZ d } 

will be called the general multiwavelet packets associated with the MRA {V 3 } of 
L 2 ( R d ) of multiplicity L 

Remark 5.1. Obviously the collection V is overcomplete in L 2 (R d ) For example 

(l) The subcollection with j = 0, n > 0, 1 5* l < L, k E Z d gives us the basic multi- 
wavelet packet basis constructed in the previous section 

(a) The subcollection with n = 1, 2, , a - 1, 1 < l < L, j E Z, k E Z d is the usual 

multiwavelet basis. 

So it will be interesting to find out other subcollections of V which form orthonormal 
bases for L 2 {R d ) Our aim is to characterize all such subcollections 
For n > 0 and j E Z, define the subspaces 

jjn = sp{ a ^ 2 ff(A 3 ■ -k) 1 <1<L, kEZ d } (5 3 °) 

Observe that . 

U°=V, and 0^ = ^. jeZ 

r=l 

Hence, the orthogonal decomposition l^+i = Vj © Wj can be written as 

^ 1 =©^ 

r=0 

We can generalize this fact to other values of u 
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Proposition 5.1. For n > 0 and j e Z, we have 


0-1 


£7 + .=®C7“ +r 


r=0 


Proof: Bv definition 


Uj + 1 = sp{a^ /, n (i4 J+1 -k) 1 < l < L, k e Z d } 


Let 


E lik {x) = a^f^A 3 * 1 -A;), for 1 < l < L, k e Z d 

Then /, A } is an oithonormal basis of the Hilbert space U™ +1 p or g < 

let 


Fl,d £ ) = E £ KmJ-MEmAx), 1 < i < L, k € 

m= 1 $££4 


and 


7f = sp{F( fc 1 < J < I, A; € 
Then, by Corollary 5 1 we have 

i?+ i=®« r 


Now 


Therefore, 


r=0 


Fl,k( X ) - £ £ Km,fi-Ak^rn,p ( x ) 

m=l ^eZ d 

m=:l a6Z d 

= E £ ^ - ■») 

m= 1 a€Z d 

= »*£ £ -*)-«) 

m= 1 a€Z d 

= a% ff n+r (A 3 x — A:), by (5 23) 


?z r = t/; n+r 


t7 + .=®^ n+r ° 

r=0 


(5 31) 


r < a — 1, 


and 
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Using Proposition 5 1 we can get various decompositions of the wavelet subspaces 
W v j > 0, which in turn will give rise to various orthonormal bases of L 2 (R d ) 

Theorem 5.2. Let j >0 Then, we have 

w > = 

= Sk* 

r=a 

(5 32) 

a ,+1 -l 

W, = © l<3 

r=a l 

w, = ® K 

r=aJ 

where U* is defined in (5 30) 

Proof: Since W 3 = U 3 , we can apply Proposition 5.1 repeatedly to get (5 32) □ 

Theoiem 5 2 can be used to construct many orthonormal bases of L 2 (R d ) We have 
the following orthogonal decomposition 

L 2 (l d ) = y 0 ®Wo©^i©^2© • 

For each j > 0, we can choose any of the decompositions of W 3 described in (5 32) For 
example, if we do not want to decompose any W 3 , then we have the usual multiwavelet 
decomposition. On the other hand, if we prefer the last decomposition in (5 32) for each 
Wj, then we get the multiwavelet packet decomposition There are other decompositions 
as well Observe that in (5 32), the lower index of U *' s are decreased by 1 in each 
succesive step If we keep some of these spaces fixed and choose to decompose others by 
using (5 31), then we get decompositions of W 3 which do not appear in (5 32) So there 

is certain interplay between the indices n € No and j € Z 

Let £ be a subset of No x Z, where No = N U {0}. Our aim is to characterize those S 
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foi which the collection 

Ps={a*/"(^ -*) 1 <1<L, keZ d , {n,j)es} 

will be an orthononnal basis of L 2 ( R?) In other words, we want to find out those subsets 
5 of N 0 x Z for which 

0 U? = L 2 (R d ) (5 33) 

By using (5 31) repeatedly, we have 

a-1 

u; = ®t/;T 

r=0 

a(n+l)-l a(n+l)~l Pa-l 

= ©£?-.= © (5 34 ) 

r=an r=an L 5=0 

a 2 (n4-l)-l a J (n+l)— 1 

= © U]- 2= ■ = © U S 

r=a 2 n r—o?n 

Let / nj = {r € N 0 o?n<r<o?{n + 1) - 1} Hence, 

Uj = © u o 

r€/n,j 

That is, 

© 17- © © ^ 

(n,j)€5 (nulls' r6/(„j) 

But we have already proved m Theorem 5.1 that 

L 2 (R d ) = © U r 0 

rG No 

Thus, for (5 33) to be true, it is necessary and sufficient that {/»,, (n,j) £ S} is a 
partition of N„ We say {A, l £ T} is a partition of * if A, C No, A,'s are pairwise 
disjoint, and U, 6 ,A, = N 0 We summarize the above discussion in the following theorem 

Theorem 5.3. let (ff . n > 0, 1 < < < L] be the bam mnttmavelet packets and 

S c N 0 x Z Then the collection of functions 

jo* ffiA* -k): 1 <1<L, keZ d , (n, j) € sj 

is an orthonormal basis of L 2 {R d ) if and only if {h,j {n,j) £ S} is a partition o/N 0 
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Let U be a separable Hilbert space A sequence {x k k € 1} of U is said to be a frame 
foi H if theie exist constants C x and C 2 , 0 < C x < C 2 < oo such that for all x € % 

CilMI 2 < E I (z.zjfe) | 2 < C 2 \\x\\ 2 (5 35) 

The largest C\ and the smallest C 2 for which (5 35) holds are called the frame bounds 
Suppose that $ = { 9 1 , <P 2 , ,<P N } C L 2 (R d ) such that {<p l ( - k) l<l<N,kel, d } 

is a frame foi its closed lineal span S(<f>). Let ip*,i/) 2 , ,ij) N be elements in £?($) So 
each i/P a huear combination of (p l ( — k), 1 < l < L, k € Z d A natural question to 
ask is When can we say that - k) 1 <1<N, fc e Z d } is also a frame for S^#)? 
If ^ € S(4>), then there exists {p 3 i k k € Z d } m l 2 ( Z d ) such that 

= E E Pjik<P l (x ~ k) 

1 = 1 kez d 

In terms of Fourier transform 

#(« = E E rne~' M m 

1=1 kez d 

= E mW(C (1<1<N), (5 36) 

1 

where P^(0 = E Pjtke~ l{k Let P(£) be the iV x iV matrix 
kez d 

Let S’ and T be two positive definite matrices of order iV We say S < T if (x, 5x) < 
(x, Tx) for all x e The following lemma is the generalization of Lemma 31m [Che] 

Lemma 5.3. Let ip l ,ip l for 1 < l < N, and, P(£) be as above Suppose that there exist 
constants C\ and C 2 , 0 < C\ < C 2 < 00 such that 

C\I < P*(Z)P(0 < C 2 I for a.e f e T* (5 37) 


Then, for all f € L 2 (R d ), we have 

CiE E K/.^t -*)>f <E E I </V( 

1 = 1 kez d 1=1 kez* 


N 


k )) | a <QE E |</V( -*)>l 

1=1 k€Z d 


(5 38) 
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Let .4 he a dilation matnx, B = A* and a = | det A\ = | det£| Let 


I<A = {a 0 ,a h ,a a _i} 


(5 39) 


and 


Kb = {P o,A, i Pa-i} (5 40 ) 

be fixed sett, of digits for A and B respectively For 0 < r, s < a - 1 and 1 < l, j < L, 
define foi ac ( 

£%{ 0 = <5(ja~ 2 e ~ l {^ +2S ~ 1 ^ 7r .“r) 


Let 


and 


E"(() = (f,7(fl) 
Btt) = 




(5 41 ) 


(5 42 ) 


(5 43 ) 

0<r,s<o-l v ' 

So E(£) is block matrix with a blocks m each row and each column, and each block is 
a squaie matnx of order L, so that E(£) is a square matrix of order aL We have the 
following lemma which will be useful for the splitting trick for frames 


Lemma 5.4. (i) If v € K A , then X) e~ t2v ( B = a£ ( 

^k b 


0v 


(n) The matnx E(£), defined tn (5 43), is unitary 

Proof: Item (z) is the orthogonal relation for the characters of the finite group Z d /BZ d 
(see [Rudl]) Observe that the mapping 

p + BZ d H- e ~ l2v { B ~ 1 i i ’ u } , ueK A 

is a character of the (finite) coset group Z d /BZ d Ifi/ = 0(ie,ifi / € AZ d ) , then 
there is nothing to prove Suppose that v ^ 0, then there exists a p € Kg such that 
g— t27r(B-W) ^ I Since k b 1S a set of digits for B, so is K B - p! Hence, 


iieK B ^ Kb 


(5 44 ) 
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Now 

£ f ,~«2jr(fl~W) _ e -t2ir (B“V',v) £ ^2* 

/'C A n h^Kb 

= jje - * by (5 44) 

msKb 

Therefore, 

£ = 0 as e -‘2*<s-VV) ± j 

To prove (n), observe that the (r,s)-th block of the matrix i? (£)£*(£) is 

E « r ‘K) (£"«))' 

£=G 

The (J,j)-th entry in this block is 

E E «&(?) (Of))' 

£-0 m=l 

= E E • ^a-VV^ 2 - 1 ^) 

t=0 m=l 

m=s l £=0 

= E Wj,A s , (by (i) of the lemma) 

m=l 

== Sij6 rs * 

This proves that E(Z)E*(£) = / Similarly, £(£)*£(£) = 1 Therefore, E(() is a unitary 
matrix D 


5.6 Splitting lemma for frame packets 


Let {<pi 1 < l < L} be functions in L 2 (R d ) such that {ipi( -k) 1 < l < L, k € Z d } 

is a frame for its closed linear span V For 0 < r < a - 1 and 1 <1 <L, suppose that 
there exist sequences {h T ljk . k € Z d } 6 l 2 (Z d ) Define f[ as m (5 4) and (5 5) That is, 

f[(x) = E E h^a^jiAx - k ) (5 45) 

]=i kei d 
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Let H t {£) be the matrix defined in (5 7) Let Ka and Kb be respectively fixed sets of 
digits for A and B as m (5 39) and (5 40) Let H(£) be the matrix 

H (0 = (ff r r(£ + 2B~ 1 /3 a n)\ (5 46) 

V ' 0<r,s<o— 1 v ' 

H(^) is a block matrix with a blocks in each row and each column, and each block is of 
order L so that //(£) is a square matrix of order aL Assume that there exist constants 
Ci and C 2 , 0 < Ci < C 2 < oo such that 


Cxi < #*(£)#(£) <C 2 I for a e £ e T* 


(5 47) 


We can wntc f{ as 


where 


//■(*) = E E h r l]k a 1/2 ipj{Ax - k) 

3= ikez d 

= E EE — a s — Ah), by (5 2) 

J=1 3=0 k£Z d 

= SEE 

J=1 i=0 k£Z d 


<p[ s \x) = a }/ 2 ipj(Ax - Q! s ), 0 < s < a - 1 


(5 48) 


Taking Fourier transform, we obtain 

L o-l 


(nrn = e t e 

3= 1 s=o kez d 

J = 1 5=0 


where pf/tf) = E 


kez d 


P"tt) = (rf/(fl) , < u < L 


(5 49) 


and 


p(e = 


(5 50) 


(5 51) 


Claim. 


H(0 = P(B()E((), 
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where E(£) is defined in (5 41)-(5 43). 

Proof of the claim The (r, s)-th block of the matrix P(Bff)E(ff) is the matrix 

e p rt mE ts (t) 

t—0 

The entry in this block is equal to 

E E 

t= 0 m=l 

= E E E 

t=o m=:i kez d 

i=o jfcgz* 

Now, the (/, j)-th entry in the (r, s)-th block of H(£) is 

h r lJ {t; + 2B- 1 M 
= a' 1 ' 2 E 

kez d 

= a-‘/ 2 E E by < 5 2 > 

= a-‘ /2 E E e " (BS, ‘ l 

1=0 ftgZ<* 

So the claim is proved In particular, we have 

= E*(Z)P*mP(BmO ( 5 52) 

Since E($ is unitary by Lemma 5 4, 2P(02f(0 P*(B0P(2tf) are similar matrices 

Let A(0 and A(£) respectively be the maximal and minimal eigenvalues of the positive 
definite matrix #*(£)#(£), and let A = inf A(fl and A = sup Afc) (It is clear from (5 51) 
that A(0 and A(£) are 27 rZ d -penodic functions ) Suppose 0 < A < A < oo Then we 
have, by (5 47) (in the sense of positive definite matrices), 

XI < H*(£)H(0 < A I for a e. £ £ T* 

which is equivalent to 

AI < P*(£)P(0 < AI forae ^ €Td 



Multiwavelet packets and frame packets of L 2 (R d ) 


79 


Then by Lemma 5 3, for all g € L 2 (R d ), we have 


a°E E E I (», ¥>!•’( -*)) 

v=0l=lk£Z d ' ' 


2 < See i (9,m -*))i 2 

s =o 1=1 kez d 

o-l L 


^ A EE E -*)) , (553) 

3=0 i=l kGZ d ' ' 


wheie <£>j s) is defined in (5 48) Since 

L : e K?.<> 1/ Vi('4 -*»i 2 =ei 

i iez* 1 3=0 i=i ftez* 1 

which follows from (5,48), inequality (5 53) can be written as 

|2 


E E \(gy'MA -*)>| 2 = EE E -*)) 2 . (sm) 

i = U 6 Z d s = 0 J = lfc e z <' 1 x ' 


AE E Ks.o'^A •-*)>!* < EE E 1(9./, '( -*)>l 2 

(=i kez d s=o z=i kez d 

< AE E |(s,« 1 'V,(A.-*)>| i 

J=i JfeeZ 1 * 


(5 55) 


This is the splitting tnck for frames the ^ -1 Z d -translates of the L dilated func- 
tions <pi(A), 1 < l < L, are ‘decomposed’ into Z d -translates of the ah functions //, 
0 < s < a - 1, 1 <1<L 

We now apply the splitting trick to the functions {// • 1 <l< L} for each s, 0 < s < 
a - 1 We have 

AEE K^W-*))! 2 * See i(9./r(-«> |J 

< ae e i<9.« i/2 /,*(At -*»i 2 . ( 5M ) 

i=i kez d 

where //’ r , 0 < r < a - 1 are defined as m (5 45) (/,' now replaces <p{) 

ff' r (x) = E E hi jk a 1/2 f;{Ax - k), 0 < s < a- 1, 1 <1 < L (5 57) 
o=i kez d 

Summing (5 56) over 0 < s < a — 1, we have 
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Using (5 55), we obtain 


A 2 £ E |(<?>a 2/ W 

1=1 kC Z d 


a-la-1 L 


- & ))l ^ E EE e \(g,f! ,r { - fc)>| 2 

»=or=oi= u ^ d '^ ,Jl v ;/l 

- l<*.a a/ W 2 -A;)) I 


(5 58) 


Now as in the case of orthonormal wavelet packets, we can define fl 1 for each n > 0 
and 1 < l < L (see (5 23) and (5 27)) In order to ensure that ff are in L 2 (R d ), it is 
sufficient to assume that all the entries in the matrix #(£), defined in (5 46), are bounded 
functions Comparing (5 57) and (5 23), we see that 


{//’ r 0 < 7 , s < a - 1} = {fi +ar *0<r,s<a — l} = 0<n<a 2 -l} 


So (5 58) can be wntten as 

a 2 E E Ks,« a/ W > -*)>| 2 < EE E l<9,/,“( -*))l 2 

1=1 kiZ d n=0 1=1 

< A 2 E E I<9,a 2/ W 2 -*)>| 2 

«=i fce z d 

By induction, we get for each j > 1 

a’E E £ e‘e E l<9,/r( -*)>l 2 

1=1 kez d n=0 1=1 fc6Zd 

< A j E £ \{g,<!?l\,(A’ -fc))| 2 (5 59) 

!=i fcezJ 

We summarize the above discussion in the following theorem. 

Note: {ff n > 0, 1 < l < L} will be called the wavelet frame packets 

Theorem 5.4. Let {ipr 1 <1<L}C L 2 (R d ) be such that {^( - k ) 1 < l < L, k e Z d } 
is a frame for its closed linear span V 0 , with frame bounds Ci and C 2 Let H(£) ,H r {£),\ 
and A be as above . Assume that all entries of H r {£ + 2B" 1 /3 s 7r) are bounded measurable 
functions such that 0 < A < A < oo. Let {f? n>0,l<l< L} be the wavelet frame 
packets and let V 3 = {/ f{A~’-) € V 0 } Then for all j > 0, the system of functions 


is a frame ofV 3 with frame bounds MCi and A J C 2 
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Proof: Since {v?i( — k) 1 < / < L, k E Z d } is a frame of Vo with frame bounds C\ and 
C 2 , it is deal that for all j 

-k) 1<1<L, keZ d } 

is a frame of \) with the same bounds So from (5 59), we have 

A- , Ci||< 7 || i < E E E \(9,fH- ~ fc ))| 2 ^ A j C 2 || 5 || 2 , for all g E V 3 (5 60) 

n=0 1=1 k£Zd 

□ 

In Thcoiem 5 1 we proved that the basic multiwavelet packets form an orthonormal 
basis for L 2 (R d ) = UV) An analogous result holds for the wavelet frame packets if the 
matrix H ( £)> defined m (5 46), is unitary 

Before proving this result let us observe how the space Uj>qVj looks like Let Vo = 
sp{(pi{ -k) 1 < l < L, k E Z d }, V 3 = {/ /(A' 1 ) € Vo} and V 3 C V 3 +i Let W = UV 3 
Then it is easy to check that f E W =£■/(• — A~ 3 k ) E W for all j E Z and k E Z d We 
claim that elements of the form A~ 3 k are dense in R d For K = {ki,k 2 , & 0 } a set of 

digits foi A, define the set 

Q =s Q(A,K) = {a: 6 x = Y, A ' 3 ^ e o e K ] 

K 3>l 

In the above representation of x, e 3 ’s need not be distinct We have 

||A-^|| < CVIMI, 

where C is a constant and 0 < a < 1 (see [Woj, Chapter 5]) For x = (xi,x 2 , , Xd) & R , 

||z|| = (Ixil 2 + |x 2 | 2 + • • + |x d | 2 )5 Therefore, the series that defines x is convergent The 

set Q satisfies the following properties (see [GM]) 


(i) Q = U? =l A~ 1 (Q + k % ) 
( u ) Ufc e z(Q + k) = R d 
(in) Q is compact 
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Let e > 0 and y € Q We first show that there exists j e 2 and k £ Z d such that 
||y - < f From (i) we have 

Q = [J A~ l {Q + k x ) 

t=i 

a r a 

=: U A (Q + k m ) + k t 

1=1 Lm=i 

= (J U + A~ 2 km + A~ l ki ) 

i=l rn=l 

Hence, foi any / > 1 and any y e Q, there exist y 3 6 Q and h,l 2 , , l 3 € K such that 

y = A~>y 3 + A~H 3 + A"**- 1 Vi + • • + A' 1 *! 

Therefore, 

||j/-4 :, {1 3 + A1 3 -i+ + A J ^ilH = ||A“' 7 y,|| 

< C'oP (as Q is compact) 

< e, choosing j suitably 

Now if y € R d , then by (ii) y = yo + p for some j/o € Q and p 6 Z d For y 0 e Q , there 
exist j > 0 and k € Z d such that ||yo - A -J fc|| < € That is, 

l|j/o + P " AT 1 (k — A J p) || < € 

=H|t,-;r ) (fc-A’p)||)||<<s 

So the claim is proved 

We have proved that W is invariant under translations by A and these elements are 
dense in R d Therefore, W is a closed translation invariant subspace of L 2 (R d ) Hence, 
W = L%{ R d ) for some E cR d (see [Rud2]), where 

L\{ R d ) = {/ e L 2 (R d ) supp f C E} 
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Now let 

L 

Eq =|J|Jb : '(supp <pi ) 

l- =lj>0 


Claim: E = Eq a e 

To piove the claim we will follow [BDR], Theorem 4 3 Since <pi(A 3 ) 6 V 3 C W, 
the function (ipi(A 3 -)) A = ± <pi(B~ 3 ) € W = {/ /6?} Therefore, B 3 (supp <pi) - 
SU pp (±ipi(B 3 )) C £ for all j > 0 and 1 < l < L, which implies that E 0 C E Let 
Ei = E \ Eq We have 




i=i 


(5 61) 


foi some 27rZ rf -periodic functions m t € L 2 ^) Hence, (5 61) implies that / = 0 on E x for 
all f eVj and hence, for all f EUV 3 = W Taking closure, we obtain that / = 0 on Ei 
for all feW But W is the set of all functions whose Fourier transform is supported in 
E Since E x C E, we get that Fi = 0ae Therefore, E = E 0 a e. □ 


Theorem 5.5. Let {^( - k) . 1 < l < L, k € Z d } C L 2 (R d ) be a frame for its closed 
linear span Vo, with frame bounds C x and C% and let Vo C V x , where V 3 — {/ f{A ) £ 
Vo} Assume that H(0 is unitary for a e £. Then {/f ( -k) n>0, l <1 < L, k eZ d } 
is a frame for the space Uj>oV) with the same frame bounds. 

More generally, let S = {{n,j) € No x be such that U (nj)es^ ls 0 P art ^ wn °f 
N 0 Then the collection of functions {a 3/2 ff(A 3 --k) \<1<L, [n,j)eS, keZ} is 

a frame for Uj>oVj with the same bounds C x and C 2 

Proof: Since H(i) is unitary, A = A = 1 so that the inequalities in (5 59) are equalities, 
and from (5 60) we have 

Cillsll 2 < "e £ E 1(9. - *=)>l 2 2 tor ° e v ‘ (5 62) 

n=0 f=l kez d 

Now let h € U 3 >oV 3 Then there exists h 3 € V 3 such that h 3 -» h as j oo Fix j, then 
for j < /, we have from (5.62) 

“] CE E I <V>/"( - *)>l 2 < c'allVll 2 

n= 0 1=1 ftez<* 



Multiwavelet packets and frame packets of L 2 ( R d ) 


84 


Letting / -Mx> first and then j oo, we have for all h € Uj> 0 V, 

EE E l(A./r( -*))l 2 <c 2 ||ft|| 2 

'*>0i=l*6Z d 

To get the i e\ erse inequality we again use (5 62) 

c.HM 2 < e'e E Ift, /,"(•-*)> I 2 

n=0 1=1 

= e'e E \(h,-h,ff( -k)) + {h,ff( -*)> I 2 

r»=0 1=1 fcgZ 11 

Theiefoie, 

C^IIMI < (E 1 EEIft-A,/r(--fe))l 2 ) i + (EEEI('>,/?'(-4))l 2 )' 

' n=0 1=1 heX* ' v n=0 1=1 keZ d ' 

< c^ii^-aii+Ce'ee i(A,/r( -*»i 2 ) i 

v n=0 1=1 k£Z d ' 

Taking j -> oo, we get 

CilHI 2 < EE E IM(-*))l a 

n^OlslfcgZ' 4 

for all /i € UVj. So the first part is proved. 

Now let Uj = sp{a 3/2 fl l (A :l - k ) : 1 <1<L, k € Z d } Then we can prove as in the 
orthogonal case (see (5 34)) that 


u? = © Vi, 

where 0 is just a direct sum not necessarily orthogonal, and I nj = {r € No o?n <r< 
a? ( n+ 1 ) _ i} Now as H(Z) is unitary, we have A = A = 1 and hence (5 56) is an equality 
Therefore, 


o—l L 


■ kez d 

Prom this we get 


E E \{s,^ff(A -*)>| 2 = EE E>r*( -*»l 

itljfeiz" r=0 1=1 keZ d 
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a 2 (n+l)-l l 

= £ ££!(«,/,'( -fc))| 2 

r=a 2 n 1=1 k€% d 

Snmlaih, 

L 0 a J (n+l)— 1 L 

£ £ /"(- 42 -* 0)| 2 = £ £ £ | ls,m - m 2 

r—o^n /= lkeh d 

= E EE 1(5, /f( ~fc))| 2 (5 63) 

r€/ nj i=l KeZ d 

Rom the first part of the theorem, we have for all / e UVj 

OII/IP < £ £ £ !</,/”( -fc))l 2 <c 2 |I/!I 2 

n> Oi=UeZ<* 

But, the set S is such that U(nj)es^ n J = Therefore, 

Cill/ll 2 < E E EE K/,/r(-fc))l 2 <^ll/ll 2 

(n j)€5 re/ n ,j «=1 kCZ d 

Using (5 63), we get 

<",u/f< £ £ £ K/.^W -*)>!’ < c» ii/ii 2 

(n,j)6S i=l AsgZ d 

for all / 6 UVj. This completes the proof of the theorem □ 
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a 2 (n+l)-l L 

= E E E 10, SI (•-*)> I 2 - 

r=a 2 n l=lk€Z d 

Similarly, 

L 0 a J (n+l)— 1 £, 

E E |<^ aI/ V”(^--(:)>| 2 = E E E l<s./r(- — *)>| 2 

1=1 k€Z d r=aJn 2 = lkeZ d 

= E EE l< 0 ,/r(--*)>l a - (5-63) 

r6/„j 2=1 fc€Z<i 

Prom the first part of the theorem, we have for all / e UVj 

CiWfW 2 <EE E \(f,n-m 2 <c 2 \\f\\ 2 . 

n>0 2=1 keZ d 

But, the set 5 is such that U (nj)es^J = ^o- Therefore, 

Ci||/ll 2 < E E E E l(/,/,’-(--*)>l 2 <C2ll/ll 2 . 

(nj)es rein, ] 2=1 kez d 

Using (5.63), we get 

Ci||/|| a < E E E \{f^ l2 f?W--k))\ 2 <c 2 \\f\\ 2 

(n,j)es 2=1 k€Z d 

for all / € UVj. This completes the proof of the theorem. □ 
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